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Abstract 



We calculate the complete next-to-next-to-leading order QCD corrections to the charm 
contribution of the rare decay K + — > -n^vv. We encounter several new features, which 
were absent in lower orders. We discuss them in detail and present the results for the two- 
loop matching conditions of the Wilson coefficients, the three-loop anomalous dimensions, 
and the two-loop matrix elements of the relevant operators that enter the next-to-next-to- 
leading order renormalization group analysis of the Z-penguin and the electroweak box con- 
tribution. The inclusion of the next-to-next-to-leading order QCD corrections leads to a sig- 
nificant reduction of the theoretical uncertainty from ±9.8% down to ±2.4% in the relevant 
parameter P C (X), implying the leftover scale uncertainties in B(K + — > -k + vv) and in the de- 
termination of \ Vtd\i sin 2/3, and 7 from the K — > nuu system to be ±1.3%, ±1.0%, ±0.006, 
and ±1.2°, respectively. For the charm quark MS mass m c (m c ) = (1.30 ± 0.05) GeV 
and \V US \ = 0.2248 the next-to-leading order value P C {X) = 0.37 ± 0.06 is modified to 
P C (X) = 0.38 ± 0.04 at the next-to-next-to-leading order level with the latter error fully 
dominated by the uncertainty in m c (m c ). We present tables for P C {X) as a function of 
m c (m c ) and a s (M z ) and a very accurate analytic formula that summarizes these two de- 
pendences as well as the dominant theoretical uncertainties. Adding the recently calculated 
long-distance contributions we find B{K + — * 7T + uu) = (8.0 ± 1.1) x 10~ n with the present 
uncertainties in m c (m c ) and the Cabibbo-Kobayashi-Maskawa elements being the domi- 
nant individual sources in the quoted error. We also emphasize that improved calculations 
of the long-distance contributions to K + — > tx + vv and of the isospin breaking corrections in 
the evaluation of the weak current matrix elements from K + —* rc Q e + v would be valuable 
in order to increase the potential of the two golden K — > -nvv decays in the search for new 
physics. 



1 Introduction 



The rare decay K + — > 7i + uu plays together with — > tt uu an outstanding role in the field 
of flavor changing neutral current (FCNC) processes both in the standard model (SM) [1] 
and in all of its extensions [2,3]. The main reason for this is its theoretical cleanness and its 
large sensitivity to short- distance QCD effects that can be systematically calculated using 
an effective field theory framework. The hadronic matrix element of this decay can be 
extracted, including isospin breaking corrections [4] , from the accurately measured leading 
semileptonic decay K + — > ir e + u, and the remaining long-distance contributions [5] turn 
out to be small [6], and in principle calculable by means of lattice QCD [7]. 

Consequently the SM decay rate of K + — > 7i + uu can be expressed almost entirely in 
terms of the Cabibbo-Kobayashi-Maskawa (CKM) [8] parameters, the top and the charm 
quark mass, and the strong coupling constant a s (M z ) that enters the QCD corrections 
calculated within renormalization group (RG) improved perturbation theory. Beyond the 
SM additional parameters like new couplings and masses of new heavy particles will be 
present in the decay rate, but from the point of view of hadronic effects, the theoretical 
cleanness of the prediction will not be affected by these non-standard contributions. 

In view of this, the theoretical uncertainties in the decay rate of K + — > tt + uu are at 
leading order essentially only of perturbative origin and in order to be able to test the 
SM and its extensions to a high degree of precision it is important to evaluate the first 
non-trivial and higher order QCD corrections to this decay mode. 

To be specific, the low-energy effective Hamiltonian for the K — > txvv system can be 
written in the SM as follows [9, 10] 

H eS = iG * a 2a Yl {KX\x c ) + X t X(x t ))(s Ll M(i7 iLl ^ L ). (1) 
V2 27r sm z v w „ 

t=e,fj,,T 

Here Gp, a, and sin 2 6^ denote the Fermi constant, the electromagnetic coupling, and the 
weak mixing angle, respectively. The sum over I extends over all lepton flavors, \ = V* s Vid 
are the relevant CKM factors and f L are left-handed fermion fields. The dependence on 
the charged lepton mass is negligible for the top quark contribution. In the charm quark 
sector this is the case only for the electron and the muon but not for the tau lepton. 

The function X(x t ) in Eq. (1) depends on the top quark MS [11] mass through x t = 
m^(fi t )/M^. It originates from Z-penguin and electroweak box diagrams with an internal 
top quark. Sample diagrams are shown in Fig. 1. As the relevant operator has a vanishing 
anomalous dimension and the energy scales involved are of the order of the electroweak 
scale or higher, the function X(x t ) can be calculated within ordinary perturbation theory. 
It is known through next-to-leading order (NLO) [10,12,13]. The inclusion of these 0(a s ) 
corrections allowed to reduce the ±6% uncertainty due to the top quark matching scale 
/it = 0(m t ) present in the leading order (LO) formula down to ±1%. Consequently the 
reached theoretical accuracy on the top quark contribution to K + — > 7r + uu and in the 
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Figure 1: Examples of Z-penguin and electroweak box diagrams that contribute both to 
K + — > tt + uu and Kl — > 7T uu. Here and in the following we do not display the neutrino 
line attached to the Z-boson. 



amplitude of Kl — ► 7r°^z/, where only X(x t ) enters, is satisfactory. 

The function X (x c ) in Eq. (1) relevant only for K + — > -n + uu depends on the charm 
quark MS mass through x c = m^(/x c )/M^. As now both high- and low-energy scales, 
namely /i w = 0(M W ) and \x c = 0{m c ) are involved, a complete RG analysis of this term is 
required. In this manner, large logarithms ln(/^///^) are resummed to all orders in a s . At 
LO such an analysis has been performed in [14]. The large scale uncertainty due to \i c of 
±26% in this result was a strong motivation for the NLO analysis of this contribution [9,10]. 

Defining the phenomenologically useful parameter 

P C (X) = ^ (^X\x c ) + V(* c )) , (2) 

with A = \V US \, one finds for A = 0.2248 at NLO [15] 1 

P C {X) = 0.369 ± 0.036 thcory ± 0.033 mc ± 0.009 aa , (3) 

where the parametric errors correspond to the ranges of the charm quark MS mass m c (m c ) 
and the strong coupling constant a s (M z ) given in Tab. 4. The theoretical error summarizes 
uncertainties due to various scales and different methods of computing a s (/x c ) from a s (M z ). 
Details on how the quoted errors have been obtained will be given in Sec. 9. 

Provided P C {X) is known with a sufficient precision, a measurement of K + — > 7i + uu, 
either alone or together with one of Kl — > tt uu, allows for precise determinations of the 
CKM parameters [1]. The comparison of this standard unitarity triangle (UT) with the one 
from 5-physics offers a stringent and unique test of the SM. In particular for the branching 
ratios B(K + — > tt + uu) and B(Kl — » n^vv) close to their SM predictions, one finds that a 

1 The numerical results presented here differ somewhat from [15]. For the present numerical evaluation 
the program used in [15] has been changed slightly in order to implement the various theoretical errors in 
an improved fashion. 
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given uncertainty a(P c (X)) translates into [15] 



\Vtd\ 
a (sin 2(3) 



sin 2(3 
(7(7) 



7 



±0.83 



±0.41 



±0.34 



(4) 



with similar formulas given in [3]. Here V t d is the element of the CKM matrix and (3 and 7 
are the angles in the standard UT. As the uncertainties in Eq. (3) coming from the charm 
quark mass and the CKM parameters should be decreased in the coming years it is also 
desirable to reduce the theoretical uncertainty in P C (X). 

To this end, we extend here the NLO analysis of P C (X) presented in [9, 10] to the next- 
to-next-to-leading order (NNLO) [15]. We encounter several new features, which were 
absent in lower orders. First, closed quark loops in gluon propagators occur, resulting in a 
novel dependence of P C {X) on the top quark mass and in non-trivial matching corrections 
at the bottom quark threshold scale fib = 0(m&). Second, the contributions from the 
vector component of the Z-boson coupling are non-trivial at NNLO and are only found 
to vanish in the sum of several contributions, which involve a flavor off-diagonal wave 
function renormalization. Third, the presence of anomalous triangle diagrams involving 
a top quark loop, two gluons, and a Z-boson makes it necessary to introduce a Chern- 
Simons operator [16, 17] in order to obtain the correct anomalous Ward identity of the 
axial- vector current [18]. The inclusion of such a Chern-Simons term is also required to 
compensate for the anomalous contributions from triangle diagrams with a bottom quark 
loop. Since all these effects arise first at NNLO, they are not included in the theoretical 
uncertainty quoted in Eq. (3), which has been estimated from the variation of scales and 
different methods of evaluating a s (fi c ) from a s (M z ). The only way to control their size is 
to compute them explicitly, which is a further strong motivation for our NNLO calculation. 

Our paper is organized as follows. In Sec. 2 we give formulas for P C {X) and B(K + — > 
at NNLO in a form suitable for phenomenological applications. In particular we 
present tables that show P C (X) for different values of a s (M z ), /i c , and m c (m c ) and we give 
a simple analytic formula for P C (X) that approximates the exact numerical result with high 
accuracy. Sec. 3 is meant to be a guide to the subsequent Sees. 4 to 7 that describe our 
calculation in detail. These sections are naturally rather technical and might be skipped 
by readers mainly interested in phenomenological applications of our result. Sec. 8 contains 
another accurate approximate formula for P C (X) that summarizes the dominant parametric 
and theoretical uncertainties. In Sec. 9 we present the numerical analysis of the NNLO 
formulas. In particular we analyze various scale uncertainties that are drastically reduced 
by going from NLO to NNLO. We present the result for B(K + — > tt + uu) and B(Kl — > ^vv) 
and we investigate the parametric and theoretical uncertainties in the determination of the 
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CKM parameters with the latter being significantly reduced through our calculation. In 
the course of this section we also present results provided by the CKMfitter Group [19] 
and the XJTfit Collaboration [20]. We conclude in Sec. 10. Some technical details as well as 
additional material has been relegated to the appendices. 



2 Master Formulas at NNLO 



2.1 Preliminaries 

In this section we will present the formula for B(K + — > 7i + uu) based on the low-energy 
effective Hamiltonian given in Eq. (1) extended to include the recently calculated contri- 
butions of dimension-eight four fermion operators generated at the charm quark scale /a c , 
and of genuine long-distance contributions which can be described within the framework 
of chiral perturbation theory [6]. These contributions can be effectively included by adding 

5P C , U = 0.04 ± 0.02 , (5) 

to the relevant parameter P C (X). The quoted error in 5P CjU can in principle be reduced by 
means of lattice QCD [7]. 

For completeness we will also recall the formula for B(Kl — > n uu) which we will use 
in Sec. 9 to obtain its updated SM value in view of the recently modified value of the top 
quark mass [21]. 



2.2 Branching Ratio for K + — > tt + w 

After summation over the three neutrino flavors the resulting branching ratio for K + 
tt + uu can be written as [6,9,10] 



B{K + -> ix + vv) 
with 



Xt ^)) + f^(x t ) + ^ (P C (X) + 5P C , U ) 



(6) 



3a 2 B(K+ -> ir°e + u) - , . u ( X V 

*+ = r K+ 1 - . ^A 8 = 5.04 ±0.17 x 10" 11 — — . (7) 

2i[ z sm v w \ 0.2248 / 

Here the parameter = 0.901±0.027 summarizes isospin breaking corrections in relating 
K + —* tt + uu to K + — > Ti°e + v [4]. The apparent strong dependence of B(K + — > tx + vv) on 
A is spurious as both P C (X) and 5P C)U are proportional to 1/A 4 . In quoting the value for 
P C (X) and B(K + -> vr+z/p) we will set A = 0.2248. a and sin 2 6 W entering B(K + -> n + uv) 
are naturally evaluated at the electroweak scale [22]. Then the leading term in the heavy 
top expansion of the electroweak two-loop corrections to X(x t ) amounts to typically — 1% 
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for the MS definition of a and sin 2 9 W [23]. In obtaining the numerical value of Eq. (7) 
we have employed a = a^{M z ) = 1/127.9, sin 2 9 W = sin 2 6 W = 0.231, and B(K + — > 
7r°e + z/) = (4.93 ±0.07) x 10~ 2 [24]. We remark that in writing B(K + -> n + vv) in the 
form of Eq. (6) we have omitted a term proportional to (X e (x c ) — X T (x c )) 2 . Its effect on 
B(K + -> ir+vis) is around +0.2%. 

The function X(x t ) entering Eqs. (1), (6) and (17) is given in NLO accuracy by 

X(x t ) = X (x t ) + ^Mx 1 (x t ) = Vx X (x t (m t )) , (8) 
4tt 

with 

r]x = 0.986 ± 0.009 . (9) 

The contribution stemming from Z-penguin and electroweak box diagrams without QCD 
corrections reads [25] 



vi \ 2 Xt + Xt 6x t -3x 2 

X {x t ) = -—^ — — ^Inxt, (10) 



X t ) 8(^1 -X t ) 



while the QCD corrections to it take the following form [10, 12, 13] 



29x t — xf — 4xf x t + 9xf — xf — xf , 

Xi{xt) = — r= Inxt 

3(1 -xt) 2 (l-xtr 



+ 8Xf + 4 ^ + "r X ' In 2 x t - 4Xt -ju 2 (l - x t ) + 8x t ^^ In ■ * 



:n) 



2(1 -x t f (l-x t y zv 11 dx t Ml 



Here Li 2 (x) = — J* dt ln(l — t)/t. The explicit /it-dependence of the last term in Eq. (11) 
cancels to the considered order the /^-dependence of the leading term X (x t (fit))- The 
factor r/x summarizes the NLO corrections. Its error has been obtained by varying fi t in 
the range 60 GeV < fi t < 240 GeV on the left-hand side of Eq. (8) while keeping /i f fixed 
at m t = m t (m t ) on the right-hand side of the same equation. The leftover /^-dependence 
in X(x t ) is slightly below ±1%. 

The uncertainty in X(x t ) is then dominated by the experimental error in the mass of 
the top quark. Converting the top quark pole mass of M t = (172.7 ± 2.9) GeV [21] at three 
loops to m t (M t ) [26] and relating m t (M t ) to mt{m t ) = (163.0 ± 2.8) GeV using one-loop 
accuracy, we find 

X(x t ) = 1.464 ±0.041. (12) 

The given uncertainty combines linearly an error of ±0.028 due to the error of m t (m t ) and 
an error of ±0.013 obtained by varying fi t in the range given above. 

As opposed to X(x t ) the charm quark contribution, represented by the parameter P C (X) 
in Eq. (2), involves several different scales like fi w , fib, and /i c . In order to control the size 
of the perturbative corrections to P C (X) the large logarithms associated with these scales 
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Figure 2: The parameter P C {X) at NNLO as a function of m c (m c ) (upper plot) and a s (M z ) 
(lower plot). In obtaining the numerical values for P C {X) we have set A = 0.2248, mt(m t ) = 
163.0 GeV, /v = 80.0 GeV, = 5.0 GeV, /i c = 1.50 GeV, a s (M z ) = 0.1187 (upper plot), 
and m c (m c ) = 1.30 GeV (lower plot). 
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P C (X) 


a s {M z ) \ /i c [GeV] 


1.0 


1.5 


2.0 


2.5 


3.0 


0.115 


0.393 


0.397 


0.395 


0.392 


0.388 


0.116 


0.389 


0.394 


0.391 


0.388 


0.383 


0.117 


0.384 


0.390 


0.387 


0.383 


0.379 


0.118 


0.380 


0.386 


0.383 


0.379 


0.374 


0.119 


0.375 


0.381 


0.379 


0.374 


0.369 


0.120 


0.370 


0.377 


0.374 


0.369 


0.364 


0.121 


0.365 


0.372 


0.369 


0.364 


0.359 


0.122 


0.359 


0.368 


0.364 


0.359 


0.354 


0.123 


0.353 


0.363 


0.359 


0.354 


0.348 



Table 1: The parameter P C {X) in NNLO approximation for various values of a s (M z ) and 
\i c . In quoting the numerical values for P C (X) we have set A = 0.2248, m t (m t ) = 163.0 GeV, 
m c {m c ) = 1.30 GeV, n w = 80.0 GeV, and fi h = 5.0 GeV. 



have to be resummed to all orders in a s using RG techniques. Keeping terms to first order 
in a s , the perturbative expansion of P C {X) has the following general structure 

p c (x) = -^pg»(x) + pP(x) + ^pwpo , (13) 

a a [fjL c ) An 

where we have suppressed the dependence of the expansion coefficients Pc k \x) on the 
involved physical and unphysical mass scales for simplicity. The leading term P^ (X) has 
been worked out in [14] while the NLO correction Pc l \x) has been calculated in [9, 10]. 

The main goal of the present paper is the calculation of the NNLO term Pc 2 \x). As 
indicated by the theoretical error in Eq. (3), the sum of the first two terms in Eq. (13) still 
exhibits sizable unphysical scale dependences, in particular the one on \i c . Besides, the NLO 
value of P C {X) depends in a non-negligible way on the method used to compute oe s (/j, c ) 
from a s (M z ) [15]. The observed numerical difference is due to higher order terms in a s and 
has to be regarded as part of the theoretical error. This source of uncertainty has not been 
taken into account in previous NLO analyses of the charm quark contribution [3,9,10]. 

As we will demonstrate in Sec. 9, the inclusion of P c (X) removes essentially the entire 
sensitivity of P C (X) on fi c and on higher order terms in a s that effect the evaluation of 
a s (fjL c ) from a s (M z ). As a result, the final theoretical error in P C (X) is reduced from ±9.8% 
at NLO down to ±2.4% at NNLO. After our calculation the theoretical accuracy on the 
charm quark contribution to K + — > tt + uu is thus also satisfactory. 
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Pc(X) 


a s (M z ) \ m c (m c ) [GeV] 


1.15 


1.20 


1.25 


1.30 


1.35 


1.40 


1.45 


0.115 


0.307 


0.336 


0.366 


0.397 


0.430 


0.463 


0.497 


0.116 


0.303 


0.332 


0.362 


0.394 


0.426 


0.459 


0.493 


0.117 


0.300 


0.329 


0.359 


0.390 


0.422 


0.455 


0.489 


0.118 


0.296 


0.325 


0.355 


0.386 


0.417 


0.450 


0.484 


0.119 


0.292 


0.321 


0.350 


0.381 


0.413 


0.446 


0.480 


0.120 


0.288 


0.316 


0.346 


0.377 


0.409 


0.441 


0.475 


0.121 


0.283 


0.312 


0.342 


0.372 


0.404 


0.437 


0.470 


0.122 


0.279 


0.307 


0.337 


0.368 


0.399 


0.432 


0.465 


0.123 


0.274 


0.303 


0.332 


0.363 


0.394 


0.426 


0.460 



Table 2: The parameter P C {X) in NNLO approximation for various values of a s (M z ) and 
iTic( m c)- 111 quoting the numerical values for P C {X) we have set A = 0.2248, m t (m t ) = 
163.0 GeV, fi w = 80.0 GeV, fi h = 5.0 GeV, and /i c = 1.50 GeV. 



The analytic formula for the sum of the first two terms in Eq. (13) can be found in [9,10]. 
The formula for P c \X) is given in Sees. 6 to 8. Setting A = 0.2248, mt{m t ) = 163.0 GeV, 
fj, w = 80.0 GeV, /ib = 5.0 GeV, and /i c = 1.50 GeV, we derive an approximate formula for 
P C {X) as a function of m c (m c ) and a s (M z ). It reads 

/ v \ 2.155 / \\ -1.417 



PAX) = 0.379 — -^rr-r y ' , (14) 

V ' \1.30GeV J V 0.H87 J K J 

and approximates the exact NNLO result with an accuracy of better than ±1.0% in the 
ranges 1.15 GeV < m c {m c ) < 1.45 GeV and 0.1150 < a s (M z ) < 0.1230. 

The dependence of P C (X) on m c (m c ) and a s (M z ) can be seen in Fig. 2, while in Tabs. 1 
and 2 we present the exact values for P C {X) for different values of a s (M z ), /i c , and m c (m c ). 
We observe that the /independence is almost negligible and the dependence on a s (M z ) is 
small. On the other hand the dependence on m c (m c ) is sizable. A reduction of the error 
in m c (m c ), which is dominated by theoretical uncertainties, is thus an important goal in 
connection with K + — > tt + uu. 

Employing the central values of the input parameters summarized in Tab. 4, the depen- 
dence of B(K + — > tx + vv) on m c (m c ) is described by 

B(A' + ^H8,5(^Y' 19 xlO- (15) 
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7 

1.2 1.3 1.4 

m c (m c ) [GeV] 




M t [GeV] 

Figure 3: The branching ratio B(K + — » 7r + z/z/) at NNLO as a function of m c (m c ) (upper 
plot) and M t (lower plot). In obtaining the numerical values for B(K + — > 7T + uu) we have 
used /i w = 80.0 GeV, //& = 5.0 GeV, and \i c = 1.50 GeV, and set all input parameters to 
their central values as given in Tab. 4. 
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which approximates the exact NNLO result with an accuracy of almost ±0.6% in the range 
1.15 GeV < m c (m c ) < 1.45 GeV. Similarly, the dependence on M t is given by 

^^ = ul (ijb) x 10 " n " (16) 

which approximates the exact NNLO result with an accuracy of better than ±0.1% in the 
range 165.0 GeV < M t < 180.0 GeV. The dependences in Eqs. (15) and (16) are exhibited 
in Fig. 3. We remark that the present analysis of the UT is practically independent of 
the exact value of the top quark mass. In obtaining both Eqs. (15) and (16) we have 
therefore set for simplicity ImA(, ReA t , and ReA c to their central values as given in Tab. 4. 
Furthermore we have used the numerical value of P C (X) evaluated at fi w = 80.0 GeV, 
fib = 5.0 GeV, and \x c = 1.5 GeV. A detailed numerical analysis of various uncertainties in 
B(K + — > it + vu) will be presented in Sec. 9. 



2.3 Branching Ratio for Kl — > 7r°vi/ 

In the case of Kl — > ti uu the charm quark contribution and the long-distance effects are 
negligible so that the relevant branching ratio is given simply as follows [9, 10] 

(TmA \ 2 

^X(x t )j , (17) 



where 



r Kh t(K l ) tn on , n nP7N w in _ 10 ( A 



k l = k + ; 7 = (2.20 ± 0.07) x 10' 



r K+ T(K+) v ' ' V - 2248 , 

Here we have summed over the three neutrino flavors and used t(Kl)/t(K + ) = 4.16±0.04 
[24]. rx L = 0.944 ± 0.028 is the isospin breaking correction from [4] with k + given in 
Eq. (7). Due to the absence of P C (X) in Eq. (17), B(Kl — > ix°vu) is plagued only by 
small theoretical uncertainties coming from fi t and Kl- The total parametric uncertainty 
stemming from mt(mt) and ImAf is on the other hand sizeable. The latter errors should be 
decreased significantly in the coming years, so that a precise prediction for B(Kl — > ifivv) 
should be possible in this decade. 



2.4 Summary 

As far as perturbative uncertainties are concerned, with the NNLO correction to the charm 
quark contribution to B(K + — > 7r + z/p) at hand, K + — > -n^vv has been put nearly on 
the same level as Kl — * ^vv. The leftover scale ambiguities are all small, so that the 
reached theoretical accuracy is now satisfactory in both decays. Similarly the errors due 
to uncertainties in a s (M z ) and m t (m t ) are small. The future of precise predictions for 
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B(K + — > tt + uu) will depend primarily on the reduction of the errors in m c (m c ) and in the 
CKM parameters, whereas B(Kl — > ir uu) is practically only affected by the uncertainties 
in the CKM elements. Non-negligible uncertainties arise in both cases also from the the- 
oretical error of the isospin breaking corrections and in the case of B(K + — > 7i + uu) from 
the uncertainty associated with the long-distance corrections. 

On the other hand the determination of the CKM parameters and of the UT from the 
K — > -kvv system will depend on the progress in the determination of m c (m c ) and the 
measurements of both branching ratios. Also a further reduction of the error in \V c b\, r^+, 
tk L i and SP CtU would be very welcome in this respect. 



3 Guide to the NNLO Calculation 

In analogy to many other FCNC processes, perturbative QCD effects lead to a sizable 
modification of the purely electroweak contribution to K + — > tx + vv by generating large 
logarithms of the form L = hi{^ 2 c / ^ v ). A suitable framework to achieve the necessary 
resummation of the logarithmic enhanced corrections is the construction of an effective 
field theory by integrating out the heavy degrees of freedom. In this context short- distance 
QCD corrections can be systematically calculated by solving the RG equation that governs 
the scale dependence of the Wilson coefficient functions of the relevant local operators built 
out of the light and massless SM fields. 

As several new features, which were absent in LO and NLO, enter the RG analysis of 
P C (X) at NNLO, the actual calculation is rather involved and it is worthwhile to outline 
first its general structure. This will be done in this section, whereas the details of the 
computation of the different ingredients that are necessary to obtain the NNLO correction 
to P C (X) will be described in Sees. 4 to 7. 

The key element of the RG analysis of the charm quark contribution to the s —>■ dvv 
transition is the mixing of the bilocal composite operators 

Q p ± = -i J AT(Q c ± (i)Q z (0) - Ql(x)Q z (0)) , (19) 

and 

Q B = -*J d 4 xT(Q c 3 (x)Ql(0) - Q u 3 (x)Ql(0)) , (20) 
where T denotes the usual time-ordering, into 

2 

Q« = -irk E (viLi^m) ■ (21) 

Here m c is the charm quark MS mass m c (/i), and the inverse powers of g have been 
introduced for later convenience, following [9]. One may arbitrarily shift such factors from 
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Figure 4: One-loop diagrams in the low-energy effective theory that involve a large leading 
logarithms. The black squares represent the insertion of the bilocal composite operators 
Q± (left) and Q B (right). The charm quark line is open in the left diagram and in all 
subsequent diagrams in which the charm quark couples to Q±. 



the Wilson coefficient into the definition of the operator. The factor fi~ 2e stems from the 
relation g = Z g gfi e , where Z g denotes the renormalization constant of g, and the fact 
that Q u written in terms of bare fields and parameters must be //-independent. All bare 
quantities will carry the subscript hereafter. The operators Q q ±, Qf,, and Q\ entering 
Eqs. (19) and (20) will be defined as we proceed. 

Interestingly the charm quark contribution to the s — > dvv amplitude involves large 
logarithms even in the absence of QCD interactions, because Q± an d Q B mix into Q v 
through one-loop diagrams containing no gluon. The relevant Feynman graphs can be 
seen in Fig. 4. Factoring out Gp and a the charm quark contribution to the s — > dvv 
amplitude thus receives corrections of 0(a™L n+1 ) at LO, of 0(a™L n ) at NLO, and of 
0(a™L n ~ l ) at NNLO. This structure of large logarithms explains the peculiar expansion 
of P C (X) in Eq. (13) with the leading term being of 0(l/a s ) rather than 0(1). 

Since there is no mixing between the bilocal composite operators Q± and Q B , the RG 
analysis of P C (X) naturally splits into two parts: one for the Z-penguin contribution which 
involves Q± and one for the electroweak box contribution that brings Q B into play. As 
the structure of the Z-penguin is more complicated than the one of the electroweak boxes 
it is useful to discuss the former type of contributions first. This will be done in Sec. 6. 
After the detailed exposition of the Z-penguin contribution it is straightforward to repeat 
the analysis in the analogous, but somewhat simpler case of the electroweak boxes. The 
essential steps of this calculation will be discussed in Sec. 7. 

The main components of the calculations of Sees. 6 and 7 performed here at the NNLO 
level are: i) the two-loop 0(a 2 ) corrections to the initial conditions of the relevant Wilson 
coefficients at /i w , ii) the three-loop 0(a>l) ADM describing the mixing of the associated 
physical operators, in) the two-loop 0{a 2 s ) threshold corrections to the Wilson coefficients 
at /if,, and iv ) the two-loop matrix elements 0(a 2 s ) of the relevant operators at /t c . 

The current-current operators Q q ± which enters the bilocal composite operator are 
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familiar from the non-leptonic | AS\ = 1 transitions. As their mixing under renormalization 
is not affected by the presence of the other operators, it is convenient to perform their RG 
analysis before discussing the Z-penguin contribution itself. This calculation involves the 
first three aforementioned steps as we will explain in Sec. 4. 

The second ingredient of the bilocal composite operator Q± is the neutral-current op- 
erator Qz, which describes the interactions of neutrinos and quarks mediated by Z-boson 
exchange. It is a linear combination of the usual vector and axial-vector couplings of the 
left-handed neutrino current to quarks, and a Chern-Simons operator that describes the 
coupling of neutrinos to two and three gluons. The inclusion of the Chern-Simons operator 
is essential to guarantee the non-renormalization of Qz to all orders in a s , which plays an 
important role in the RG analysis of the Z-penguin sector. The subleties arising in con- 
nection with Qz will be reviewed in Sec. 5.1 before analyzing the Z-penguin contribution 
itself. 

Finally, we will discuss the operators Q\ and Q\. These operators are the building 
blocks of the bilocal composite operator Q B and describe the interactions between leptons 
and quarks mediated by W^-exchange. They will be briefly discussed in Sec. 5.2. 

In summary, after the three preparatory sections, namely Sees. 4, 5.1, and 5.2, that 
discuss the dimension-six operators Q q ±, Qz, Ql, and Q% the actual NNLO calculation 
relevant for the charm quark contribution to K + — > 7r + uu is presented in Sees. 6 and 7 
for the Z-penguin and the electroweak box contributions, respectively. The result of these 
efforts will be summarized in Sec. 8. 



4 Current-Current Interactions 
4.1 Effective Hamiltonian 

As we are only interested in the charm quark contribution to the s — > dvv transition, 
we can drop the parts of the low-energy effective Hamiltonian that are proportional to 
the CKM factor A t . The unitarity of the CKM matrix then allows one to express all 
the relevant contributions in terms of one independent CKM factor, namely A c . P C (X) 
receives contributions from Z-penguin and electroweak box diagrams with internal charm 
and up quarks. Examples are depicted in Fig. 1. For scales /i in the range fi c < /i < fi w 
the four-quark interaction mediated by iy-boson exchange is described by the effective 
current-current Hamiltonian 

AC 

= -^A c £ C 4 (ai) {Q1 - Qt) , (22) 

» i=± 

where 

Ql = \ {{stlMilO ± {sll^Mrdt)) . (23) 
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Figure 5: Sample diagrams for the 0{a 2 s ) corrections to the initial values of the Wilson 
coefficients of the current-current operators. 

Here a and j3 are color indices. At LO the operators in Eq. (23) renormalize multiplicatively. 
Beyond LO they mix into so-called evanescent operators, which vanish algebraically in 
n = 4 dimensions [27-30], but affect the values of the Wilson coefficients C±(fi). These 
operators can be chosen in such a way that the renormalized matrix elements of Q q ± and 
its Fierz transform are the same. For this choice Q q ± have well-defined and distinct isospin 
quantum numbers and do not mix with each other at NLO and beyond. The definitions of 
the evanescent operators required to preserve the diagonal form of the NNLO anomalous 
dimension matrix (ADM) in the Q± basis can be found in App. A.l. 

4.2 Initial Conditions 

We now turn our attention to the calculation of the initial conditions of Q q ±. Dropping the 
unnecessary flavor index q we write for i = ± 

CM = cf'Ov) + ^C®M + (^) 2 Cf\» w ) , (24) 

where a s (/j, w ) denotes the strong coupling constant in the MS scheme for five active quark 
flavors. The values of the coefficients C± (;Uw) are determined by matching Green's func- 
tions in the full and the effective theory at \i w . In the NNLO approximation this requires 
the calculation of one-particle-irreducible two-loop diagrams. Sample SM graphs are dis- 
played in Fig. 5. For what concerns the regularization of infrared (IR) divergences we 
follow the procedure outlined for example in [31,32], which consists in using dimensional 
regularization for both IR and ultraviolet (UV) divergences. While the former singularities 
are removed by renormalization, the latter poles cancel out in the difference between the 
full and the effective theory amplitudes. For detailed descriptions of higher-order matching 
calculations of strong and electroweak corrections applying the latter method we refer the 
interested reader to [32-34]. 

Using naive dimensional regularization (NDR) [35] with a fully anticommuting 75, we 
obtain for the standard choices of the Casimir invariants Ca = 3, Cp = 4/3, and five active 
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quark flavors 



21 '3 Mi 



1 . o 1 / 1 



4 2) (/v) = - ggQo (135677 =F 124095) + - (7 ± 51) tt 2 T - ( 1 T - ) 



(25) 



4(HT249)ln^ + i(7±51)ln 2 ■ 



w 



36 v ' M£ 6 V 

The function T(x t ) depends on the top quark MS mass via x t = m%((J, w ) / M%y- It originates 
from SM diagrams like the one shown on the left of Fig. 5. Subtracting the corresponding 
terms in the gluon propagator in the momentum-space subtraction scheme at q 2 = 0, which 
ensures that a s (fi w ) is continuous at the top quark threshold fi w , 2 we find 



112 /20 
T{x t ) = — + 32x t + f — + 16x t ) lnx t 



(8 + 16x t ) V^x t - 1 Cl 2 ^2 



arcsin 



(26) 



where Cl2(x) = Im(Li2(e* x )). As far as the one-loop initial conditions, namely C± (l^w) are 
concerned, our results agree with those of [27,36]. They also agree with the results obtained 
in [37,38] after a transformation to our renormalization scheme specified by Eqs. (A. 2). The 
general formalism of a change of renormalization scheme discussed in detail in [39] can also 
be used to verify that our result for the two-loop initial conditions C± (/%-) coincides with 
the findings of [32]. This is shown in App. A.2. 



4.3 Anomalous Dimensions 

The Wilson coefficients are evolved from fi w down to the relevant low-energy scale n with 
the help of the RG equation. In this way, large logarithms of the form ln(/i 2 / '/i^-) are 
resummed to all orders in a s . In mass-independent renormalization schemes like MS the 
RG equation is given by 

H^-CiO*) = yMCjQi) , (27) 

where Jji(^) is the entry of the ADM describing the mixing of Qj into Qi. In the case of 
Q q ± we will denote the diagonal entries of 7(/i) by 7±(/^). 

2 This scheme coincides with MS for fi < fx w . For details see for example [32,34]. 
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Figure 6: Sample diagrams for the 0(a^) mixing among the current-current operators. 
The fermion loop in the right Feynman graph contains / active quark flavors. 



In the NNLO approximation the ADM has the following perturbative expansion 



"s(X> ~(o) , / a s(v) \ 2 , ( a s(v>) \ 3 £ 



w = i^H*J + ' + \^) f > (28) 

where the coefficients 7^ can be extracted from the one-, two-, and three-loop QCD 
renormalization constants in the effective theory. The renormalization matrices are found 
by calculating amputated Green's functions with single insertions of Q q ± up to three loops. 
Sample diagrams are shown in Fig. 6. The corresponding amplitudes are evaluated using 
the method that has been described in [40-42]. We keep the gauge parameter arbitrary 
and find it to cancel from 7^, which provides a powerful check of our calculation. To 
distinguish between IR and UV divergences, we introduce a common mass M for all fields 
and expand all loop integrals in inverse powers of M. This makes the calculation of the 
UV divergences possible at three loops, as M becomes the only relevant internal scale, and 
three-loop tadpole integrals with a single non-zero mass are known [41,43]. Comprehensive 
discussions of the technical details of the renormalization of the effective theory and the 
actual calculation of the operator mixing are given in [40,42]. 

While 7± is renormalization-scheme independent, 7^ and 7^ are not. In the NDR 
scheme supplemented by the definition of evanescent operators given in Eqs. (A. 2), we 
obtain for Ca = 3, Cp = 4/3, and an arbitrary number of active quark flavors / 



7i° } = ±6 I 1 T 3 
7i 2) = (349049 ± 201485) — (115577 q= 9795) / 



1 
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Here ((x) is the Riemann zeta function with the value £(3) ~ 1.20206. Again, we find 
agreement with the one- and two- loop results of [27,36], and therefore also with the findings 
of [37,38,42] that were obtained in different renormalization schemes. We also confirm the 
three-loop results presented recently [39]. The explicit formulas that allow the conversion 
of the latter anomalous dimensions to our scheme can be found in App. A.2. 



4.4 Threshold Corrections 

In order to compute the Wilson coefficients for scales \i much lower than fi w a proper 
matching between effective theories containing / and / — 1 active quark flavors has to be 
performed each time one passes through a flavor threshold. It is achieved by requiring 
that the Green's functions in both effective theories are the same at the point [if = 0(rrif) 
where the quark with mass m/ is integrated out. This equality translates into 

{Q f '\^)) T Cf-\^) = {Q f {n f )f&{ N ) , (30) 

where the labels / and / — 1 indicate to which effective theory the variable belongs. 
Including corrections up to NNLO we write 

(Q'Mf = <<5'> ,0,T (i + + 2 f ■ (3D 

where (Q^)^°\ r^ 1 -* (///), and r^ 2 \[if) codify the tree-level, one-, and two-loop matrix 
element of the column vector Q* containing the relevant physical operators. The other 
quantities entering Eq. (30) can be expanded in a similar fashion. 

Another subtlety arises when working in mass-independent renormalization schemes, 
because the matching conditions connecting the strong coupling constants of the effective 
theories with / and / — 1 active quark flavors are non-trivial in such schemes. In particular, 
in the MS scheme one has in the NNLO approximation [44, 45] 3 




An 3 mj 



(32) 



where rrif = m/(m/) denotes the quark MS mass. Of course, the appearance of both 
logarithmic and finite corrections can be avoided by adjusting the renormalization scheme 




3 It should be noted that the non-logarithmic term in Eq. (32) differs from the result published in [44]. 
The authors of [44] have revised their original analysis and have found agreement with [45]. 
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Figure 7: Examples of diagrams that contribute to the 0(a s ) matrix elements of the 
current-current operators. 



so that ct{ _1 (/iy) = a{ (///). While this does not affect the physical amplitudes [39], one 
leaves the class of mass-independent renormalization schemes with the drawback that the 
usual RG equations do not hold below the matching point and the resummation of large 
logarithms gets obscured. Hence it is much more convenient to stick to the MS prescription 
of a s and to apply Eq. (32) whenever a flavor threshold is crossed. We will follow this 
approach below. 

In terms of the discontinuities 

sd^M = c^ k \ N ) - df-Mfrf) , #<*>(/*/) = fw>(/*/) - tf-WM , (33) 

the solution of Eq. (30) can be written in a relatively compact form. Up to the second 
power in the strong coupling constant we obtain 

*# 0) (A*/) = 0, 



= - (srWfrf) + |ln^ (34) 

where the second line resembles the NLO result derived in [46]. Note that at NNLO the 
logarithmic 0(a 2 s ) correction entering the right-hand side of Eq. (32) starts to contribute 
to the matching conditions of the Wilson coefficients at each flavor threshold. 

In the case of C±(p) the explicit expressions for the threshold corrections turn out to 
be much simpler than suggested by Eqs. (34), because the corrections 5r±\nf) vanish, as 
the 0(a s ) matrix elements of the current- current operators are identical in the effective 
theories. Sample diagrams can be seen in Fig. 7. In consequence 5C±\nf) = 0. Note that 
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Figure 8: Sample diagrams for the 0(a 2 s ) matching corrections to the Wilson coefficients 
of the current- current operators at each flavor threshold. At these thresholds the number 
of active quarks changes from / to / — 1. 



this is in contrast to the case of the QCD and electroweak penguin operators which receive 
non-trivial threshold corrections at NLO [46,47]. However, non- vanishing discontinuities 
8r±\/j,f) arise from the diagrams depicted in Fig. 8. For what concerns the calculation of 
the graphs itself, we have adopted two different methods to regulate IR divergences and 
found identical results for the threshold corrections. The first approach mentioned earlier, 
uses dimensional regularization for both IR and UV singularities and calculates on-shell 
matrix elements with zero external momenta. The second method uses small quark masses 
as IR regulators and computes matrix elements with zero external momenta which are 
now off-shell. Useful details on the latter procedure can be found in [48]. The unphysical 
coefficients f^ k \fif) differ in both cases and depend on the IR regulators. However, this 

(2) 

dependence cancels in the combination entering 5C± (///). The correct implementation of 
the discontinuity in a s of Eq. (32) and of similar decoupling relations for the gluon and 
quark fields are of crucial importance for this cancellation. In the second method using 
off-shell matrix elements another subtlety occurs, as now both sides of Eq. (30) depend on 
the gauge parameter, and in order to obtain a gauge-independent and IR-safe result one 
has to take into account that the gauge parameter is discontinous across flavor thresholds 
as well. We do not give the decoupling relations for the gluon, the quark field, and the 
gauge parameter here. They can be found for example in [45]. 

At the bottom quark threshold scale /j& we find for the non-trivial matching conditions 
of the Wilson coefficients of the current- current operators in the NDR scheme 

*H)(I+H+^£S)). 

where r\ h = a s (fi w ) / a s (fib) and m b = m&(m&) denotes the bottom quark MS mass. We also 
remark that diagrams like the one shown on the right of Figs. 7 and 8, which correspond to 
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QCD corrections to a current operator, do not contribute once all counterterms have been 
included. The coefficients 5C±\fib) depend on the renormalization scheme chosen for Q±, 
in particular on the specific structure of the evanescent operator defined in the second line 
of Eqs. (A. 2). Choices other than these would lead to operator mixing between Q q ± . 

The discontinuities 5C±\/j c ) at the charm quark threshold scale /i c can be ignored as 
it is more convenient to express the final low-energy Wilson coefficient in terms of the a s 
of the effective theory with four active quark flavors rather than in terms of the a s of the 
effective theory with three active quark flavors, because no RG equation needs to be solved 
below /i c . This will be explained in more detail at the end of Sec. 6. 



5 Neutral and Charged Currents 
5.1 Neutral Current: Z-boson Exchange 

The low-energy effective Hamiltonian describing the interactions of neutrinos and quarks 
mediated by Z-boson exchange is given by 

W -=M^ CzM0 - (36) 

where 

Qz = E ( ( J ? - 2e « sin2 M ®v - $ + Qcs)), (37) 

and the sum over q extends over all active light quark flavors at the renormalization scale 
/i, while Ig = (+1/2,-1/2) and e q = (+2/3,-1/3) denote the third component of the 
weak-isospin and the electric charge of the up- and down-type quarks, respectively. The 
appropriate normalization of the electromagnetic coupling a and the weak mixing angle 
sin 2 6 W will become clear after our discussion in Sec. 6. 

Removing the Z-boson as an active degree of freedom from the effective theory induces 
a vector as well as an axial-vector coupling of the left-handed neutrino current to quarks 

Qv= E (riMviiffvii) , Qa= E (<?7/.75<?)(^l7 / ^l) • (38) 

£=e,fi,r l=e,fx,T 

In the literature on K + — > 7r + z/P the operator Qy is usually omitted, 4 as is does not 
contribute to the decay rate through NLO. We keep Q v throughout our NNLO calculation. 
While individual diagrams are non-vanishing, we verify explicitly that both the two-loop 
matching diagrams and the three-loop mixing diagrams with Q v sum to zero. We will 
discuss this issue in more detail in Sec. 6 after presenting our final results for the anomalous 
dimensions and matrix elements, respectively. 

4 To our knowledge the only exception is the recent publication [6]. 
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Figure 9: The one- loop top quark contribution to the coupling of the Z-boson to two and 
three gluons. Diagrams obtained by the interchange of the external gluons are not shown. 



As can be inferred from Fig. 9, decoupling the top quark generates furthermore an 
effective gauge- variant coupling of the left-handed neutrino current to two and three gluons 
which can be expressed in terms of the following Chern-Simons operator 

Qcs = £~ 2 e™" 4 (G a ^d, 2 G; 3 + lgf abc G^ 2 G^ £ {y lLl ^ L ) . (39) 

Here g denotes the strong coupling constant, e^ 1 ^ 2 ^ 3 ^ is the fully antisymmetric Levi-Civita 
tensor defined with e 0123 = +1, is the gluon field, and f abc are the totally antisymmetric 
structure constants of SU(3). We remark that the 't Hooft-Veltman (HV) prescription [49] 
and dimensional reduction (DRED) [50] lead to the same result for the triangle diagrams, if 
the mathematically consistent formulation of the DRED scheme presented recently in [51] 
is employed. A description of the HV scheme can be found for example in [27]. 

The inclusion of Qcs in Eq. (37) is essential to obtain the correct anomalous Ward 
identity for the axial- vector current [18] and in consequence to guarantee the vanishing of 
the anomalous dimension of Qz to all orders in perturbation theory. We stress that we do 
not add this contribution in an ad hoc way, instead Qcs is generated in an unambigous 
way from the diagrams in Fig. 9. Our effective theory is anomaly free, because Qcs cancels 
the anomalous contribution from the triangle graph with a bottom quark, just as the 
anomalous effects from top and bottom quarks cancel in the SM. Let us illustrate how 
the cancellation between the contributions from Q\ and Qcs to the anomalous dimension 
of Qz occurs at lowest order. As depicted in Fig. 10, the first non-trivial mixing arises 
at 0(a 2 s ) with Q q A mixing into itself through two-loop diagrams and Qcs mixing into Q\ 
through a one-loop diagram. Choosing the operator basis as (Qz,Qcs), we find for the 
NLO anomalous dimension matrix 5 

f = Cf (-12 3 • (40) 

5 The calculation has been performed in the background field gauge for the gluon field [52], which makes 
it possible maintain explicit gauge invariance at the level of off-shell Green's functions [30], keeping the 
gauge parameter arbitrary. 
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in agreement with [17]. At NNLO one needs further an evanescent operator so that 7« is 
enlarged to a 3 x 3 matrix. 

In the chosen operator basis the LO contributions to the initial values of the Wilson 
coefficients are 

(CzM,C C sM) = (1,0). (41) 

The particular form of then ensures that the Wilson coefficients up to NLO satisfy 
{Cz(ij), CcsilA) = (1) 0) a t an y scale \i. In fact, this scale independence is a striking conse- 
quence of the Adler-Bardeen theorem [53], which states that the Adler-Bell-Jackiw (ABJ) 
anomaly [18] of the axial- vector current is not renormalized in perturbation theory. This 
theorem is strictly proven to all orders for the abelian case [30,53], while strong arguments 
suggest that it holds true for the non-abelian case [54]. Assuming that the ABJ anomaly 
equation survives renormalization, it is easy to show that Cz{fj) is scale independent if 
and only if Ccs{v>) does not receive radiative corrections in the chosen operator basis, 
where Eq. (41) holds. In a renormalizable anomaly-free theory, such as the SM [55], this 
can always be achieved by invoking an additional finite renormalization of the axial-vector 
current [17,30,56]. For what concerns Czip) this means that one has to perform a fi- 
nite renormalization of Q\ to obtain the matching condition Ccs(^w) — beyond one 
loop. The corresponding finite 0(a 2 s ) correction to the renormalization constant Zaa will 
be computed in Sec. 6. Also this finite renormalization is not an ad hoc addition to our 
calculation, but originates from the loop diagrams in Fig. 17 containing a top quark. In- 
stead of absorbing these effects into Zaa one could include them in Cz(^w)- in this case 
one would also find a non-zero anomalous dimension of Qz- Since both terms combine 
to reproduce the effect of Zaa the physical result is however unchanged. Furthermore, 
owing to the definition of Qz in Eq. (36), the Wilson coefficient of Qcs does not receive a 
matching correction at the bottom quark threshold scale fib- As the anomalous dimensions 
of operators in the effective theory correspond to coefficients of large logarithms in the 
full theory, the RG invariance of Cz{fi) implies that anomalous subdiagrams involving the 
Z-boson do not give rise to logarithms ln(m^/mf) in the associated SM amplitudes to all 
orders in perturbation theory. In contrast, large logarithms proportional to m^/M^, which 
are relevant to our calculation, may arise. Such terms correspond to higher- dimensional 
operators, which are a priori not covered by the Adler-Bardeen theorem. In Sec. 6 we will, 
however, show by an explicit three-loop calculation that the latter terms are absent in the 
0(a 2 s ) charm quark contribution to the s — > dvv transition in the SM. 

5.2 Charged Current: W- Exchange 

In contrast to the neutral-current case, the discussion of the effective charged-current cou- 
plings can be kept rather short. The interactions between leptons and quarks mediated by 
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Figure 10: Diagrams describing the mixing of Q\ and Qcs m t° Q\ an d E q A at 0(a 2 s ). The 
fermion loop containing / quark flavors is anomalous, if / is odd. The Feynman graph in 
which the fermion flow in the closed quark line is opposite to the one shown in the left 
diagram is not displayed. 



Vy-boson exchange are encoded in the following low-energy effective Hamiltonian 

AC 1 

H% = ~7§ Cw{») £ (V qs Ql + V q * d Ql) , (42) 

q=u,c 

where 

Ql= £ (^7^)(^7 M 4), Ql= Yl {wMtfrfvii.) ■ (43) 

£=e,fj,,T £=e,fi,T 

Since the effective charged-current couplings Q q z and Q\ do not mix under renormalization 
the Wilson coefficient CV(/i) is /i- independent. The normalization of Eq. (42) is chosen 
such that Cw{p) = 1- 



6 Z^-Penguin Contributions 
6.1 Effective Hamiltonian 

After integrating out the top quark and the heavy electroweak gauge bosons we first en- 
counter an effective Hamiltonian which is valid for scales \i in the range \i c < fi < \x w with 
dynamical bottom and charm quark fields. The Z-penguin contribution involves Ti,^ and 
7i^ s defined in Eqs. (22) and (36) as well as the effective Hamiltonian 

K P ff = %A c TO C?(jj.)Q v , (44) 
V 2 sm z 6 W 

which finally brings the leading dimension-eight operator Q v of Eq. (21) into play. 

The desired matrix element (tt + uu\T\K + ) involves the transition operator T = T p + 
T B , where T p and T B denote the Z-penguin and electroweak box contribution, respec- 
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Figure 11: Examples of Z-penguin diagrams that contribute to the initial value of the 
Wilson coefficient of the leading dimension-eight operator Q v at 0(a 2 ). 



tively. The Z-penguin contribution to the transition operator takes the following form 



Notice that in passing from the first to the second line we have used Cz{^) = 1- The last 
two terms in Eqs. (45) are the bilocal composite operators Q± that involve the effective 
current- and neutral-current couplings Q± and Qz- The former operator has already been 
introduced in Eq. (19). In contrast to [9] we have defined it in terms of chiral and not 
"V — A" fermion fields. This results in the factors of 4 multiplying C±(/i) in the above 
equation. 

As the normalization of a is determined by the short distance interactions at fi w , 
the Wilson coefficients are appropriately expressed in terms of Gf using the relation 
g 2 (fi w )/4:TT = a(/J, w )/ sin 2 9 w (fj, w ) = y/2/ir GpM^-, where all running couplings are de- 
fined in the MS scheme. The typical case is that of electroweak box diagrams, to be 
discussed in the next section, which are clearly proportional to g 4 (/x w )/M^ oc G 2 F M^. 
After decoupling, all short distance information is encoded in the Wilson coefficients and 
in Gf, which does not evolve in the effective theory. Hence the electromagnetic coupling 
a and the weak mixing angle sin 2 9 W entering Eqs. (45) and (93) are naturally evaluated 
at the weak scale [22]. 

6.2 Initial Conditions 

The matching corrections for (/i) are again found by requiring equality of perturbative 
amplitudes generated by the full and the effective theory. Examples of two-loop Z-penguin 
SM diagrams can be seen in Fig. 11. Regulating spurious IR divergences dimensionally we 




(45) 



Gf 7t« 



(Cfo)Q„ + 4C + (j*)Qj + 4C-(n)Q*) . 



v/2 c Ml sin 2 6 



w 
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Figure 12: Diagrams involving the double operator insertions (Q±,Q q A ), (Q±,E\), and 

(QlQD- 



obtain for the non-zero matching conditions in the NDR scheme 



(46) 



(fi w ) = 4CV 33 + 4vr 2 + 34 In ^- + 12 In 



2 




w 



where the first line recalls the NLO result [9], while the second one represents the new 
NNLO expression. 

6.3 Anomalous Dimensions: Non- Anomalous Contributions 

The mixing of lower- into higher- dimensional operators through double insertions leads in 
general to inhomogeneous RG equations [29]. In the case of the Wilson coefficient C^(fi) 
introduced in Eq. (45) one has explicitly 



where 7„(/i) encodes the self-mixing of Q u , while the anomalous dimension tensor 7^ v {n) 
describes the mixing of the bilocal composite operators Q± into Q u . The factor of 4 in the 
above relation is of course a direct consequence of the factors of 4 in Eq. (45). 

Since the conserved current Sl7/^l in Q u is not renormalized, the RG evolution of Q u 
stems solely from the prefactor rn^/g 2 in Eq. (21). In terms of the expansion coefficients 
of the anomalous dimension of the charm quark MS mass and of the QCD /3-function, the 
corresponding anomalous dimension reads 



(47) 



7, 



'?> = 2 (7<? " Pk) ■ 



(48) 
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Figure 13: Sample diagrams for the mixing of the double insertion {Q±,Qz) m t° Qu & t 



In the particular case of QCD one has up to the NNLO level 6 



7 <o>.8, = = 2498 - (If + f><(3)) / - f?/ 



a 11 2 f a mo 38 / a 2857 5033 , , 325 f 2 



(49) 



The contributions to the anomalous dimension tensor 7±„(ju) stemming from non- 
anomalous diagrams can be decomposed in the following way 



7 ± T = -i^^,^ , -(i-| S m 2 » l ,)^S , 



(50) 



where the superscript A and V marks the corrections arising from diagrams with a double 
operator insertion (Q±,Q\) and (Q±,Q V ). In the NDR scheme supplemented by the 
definition of evanescent operators given in Eqs. (A. 2) we obtain after setting Ca = 3 and 
Cp = 4/3 the following coefficients 

7£S? = -4 (1 ± 3) , 

7^? = 16 (2 7±S?=0, (51) 

lif = (45124 ± 484917) + 32 (13 ± 15) C(3) ± 144/ , 

The results in the second line agree with the findings for 7^3 of the prior NLO calculation [9] 

if one takes into account i) a factor of —1/2 arising from the decomposition of 7±^ in 
Eq. (50) and ii) a factor of 4 that stems from the fact that our operators Q± and Q v are 

6 We have calculated the anomalous dimension of the quark mass and the strong coupling constant in 
the MS scheme up to the three-loop level, finding perfect agreement with the literature [17,57]. 



V(0) 

i±y 


= 0, 


V(l) 

i±y 


= 0, 


V(2) 

7±,„ 


= 0. 
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Figure 14: Sample diagrams for the O(a^) mixing of (Q±,Qz) into Q 



defined in terms of chiral and not as traditionally done U V— A" fermion fields. The third line 
shows our new NNLO results. We stress that also at NNLO the part of the double operator 
insertion (Q±,Qz) proportional to Q±Q C A — Q±Q U A accounts for the complete mixing. In 
order to understand this feature it is important to realize that one can distinguish two 
kinds of contributions: i) diagrams in which Qz couples to an up-type quark as on the 
left of Figs. 12, 13, and 14, and ii) diagrams in which Qz couples to a down-type quark 
as on the right of Figs. 13 and 14. This classification holds true to all orders in QCD. 
Diagrams of type i) containing an insertion of the vector part of Qz do not contribute 
to the anomalous dimensions. Further diagrams of type ii) vanish in LO and are UV- 
finite at NLO, but do have UV poles at NNLO. However, their overall contribution is 
cancelled by diagrams like the ones shown in Fig. 15, which induce a flavor off-diagonal wave 
function renormalization. Of course, this additional wave function renormalization has to 
be included in the renormalization of Qz- Finally let us mention that the contributions 
from diagrams of type ii) containing an insertion of the axial- and vector part of Qz have 
opposite signs, which is a consequence of Q q ± containing only left-handed down and strange 
quark fields. 

6.4 Anomalous Dimensions: Anomalous Contributions 

In order to calculate the 0(ce^) mixing of the double insertions (Q±, Qz) info Q u , we will 
need the renormalization constants of Q q A and Qcs defined in Eqs. (38) and (39) up to 
0(a%), because these operators appear as subloop counterterms in the effective theory. 
Diagrams involving such subgraphs can be seen in Fig. 16. The renormalized operators Q\ 
and Qcs defined in Eqs. (38) and (39) may be expressed in terms of the bare ones as 

Qa = ZaaQa o + ZaeEa o > 

-2e q , ( 52 ) 

Qcs = A 4 e Qcs,o + Z C aQa,o + %ceE A q , 

where the unexpected factor /i _2e stems from the relation between the bare and the renor- 
malized strong coupling constant. Note first that Q q A is protected by its gauge invariance 
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s q c d d d s Q C d d d 



Figure 15: Sample diagrams for the flavor off-diagonal wave function renormalization at 

from non-diagonal renormalization involving the gauge- variant operator Qcs, an d second 
that Qcs has no diagonal renormalization due to the factor g 2 in its definition. 

Within dimensional regularization the renormalization of Q q A is not exhausted by a 
multiplicative factor, but involves mixing with the following evanescent operator as well 

& A = i^^(rtpftnW) E feWtt) + «%- (53) 

£=e,/j.,T 

In our case the quark fields in Q\ and E q A correspond to an open fermion line. We will 
calculate the parts of the anomalous dimensions that involve anomalous subloops with 
insertions of Q q A using three different prescriptions for 75 for this open line, namely NDR, 
HV, and DRED. Together with the two possibilities to treat 75 in the closed fermion 
loop, which are HV and DRED, this amounts to six renormalization prescriptions in total. 
It is instructive to see how different scheme-dependent pieces combine into a scheme- 
independent result for P C (X). Note that diagrams with an insertion E q A must be included 
not only in the NDR and HV schemes [27-30], but also in DRED, which is unexpected at 
first sight. The crucial point here is that a mathematically consistent definition of DRED 
involves infinite-dimensional spaces just as NDR and HV: the DRED scheme entails i) a 
formally 4-dimensional, but really infinite-dimensional, space for the gauge fields and Dirac 
matrices, and ii) a formally n-dimensional space for the momenta, which is a subspace of 
the former one [51]. In consequence, E q A is not identical to zero in DRED, as it belongs to 
the formally n — 4-dimensional complement of the n-dimensional space. 

The renormalization constants Z^a^ Zae, Zca-, and Zce entering Eqs. (52) are found by 
calculating the UV divergent parts of Feynman diagrams in the effective theory. Sample 
graphs encoding the first non-trivial mixing of the set (Q% Qcs) into Q\ and E A are 
displayed in Fig. 10. In the MS scheme we obtain 

7(2,1) _ or- 7(2,1) _ I r 

^aa — oof, Zj ae ——-^f, 

z (54) 

7(2,1) 7(2,1) 

Zj ca ——ou f , Zj ce —u f , 

where the symbol Z^'^ denotes the coefficient of the l/e l pole of the 0(a k s ) term of the 
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Figure 16: Examples of diagrams involving the double operator insertion (Q±,Q q A ) and 
(Q±, Qcs)- in the (Q±, Q 9 A ) case the fermion loop containing / quark flavors is anomalous, 
if / is odd. 



associated renormalization constant. By taking into account the factor fi~ 2e present in 
the second line of Eqs. (52) and switching to the basis (Qz,Qcs) one recovers the ADM 
given in Eq. (40). All the one-, two-, and three-loop results presented in Eqs. (54), (55), 
(62) and (63) are again calculated using two different methods. In the first approach, IR 
singularities are regulated by introducing a common mass parameter into all the propagator 
denominators including the gluon ones [40,41], while in the second one only the mass of 
the open quark line is kept non-zero [48]. The two methods give the same results for the 
MS renormalization constants. 

The aforementioned finite renormalization of Q\ is most easily found by insisting that 
the 0{a 2 s ) correction to its initial condition is identical to zero. This matching requires the 
calculation of the two graphs shown in Fig. 17. The UV divergences from these diagrams 
are canceled by a counterterm proportional to Z AA \ The leftover finite contribution can 
be either absorbed into the initial value of the Wilson coefficient of Q q A or into a finite 
renormalization constant Z AA \ The latter possibility is more convenient, as it avoids a 
spurious RG running of the effective neutral-current coupling Qz, which would otherwise 
occur beyond NLO. We find 

Z^ = \C F , (55) 

for both HV and DRED. In the former scheme we reproduce the result of [17]. With 
the superscript A we indicate that we have only considered the contributions related to 
anomalous graphs. The contributions from all other diagrams, which have only open 
fermion lines, can be discussed separately. If the HV scheme is used for the latter diagrams, 
an additional finite renormalization constant is needed [17,30,56]. 

While we could rely on the Adler-Bardeen theorem to ensure that the contributions 
of Q\ and Qcs cancel in the renormalization of Qz, the situation is more complicated 
in the case of the transition operator encountered in Eq. (45), because it involves double 
insertions of the type (<3±,<3^) and (Q±, Qcs)- Typical examples of such graphs are 
shown in Fig. 16. Apparently, their renormalization requires counterterms proportional to 
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Figure 17: The SM diagrams one has to compute in order to find the 0{a 2 s ) correction to 
the initial value of the Wilson coefficient of Q q A . 



Q u . The associated renormalization constants Z± • (jii) can be extracted at any given order 
in a s by requiring 

(Q^W + z^MiQM), (56) 

to be UV finite. Here 

Q±i = ~ij d*xT{Q c ± {x)Q 3 {0) - Ql(x)QM) , (57) 

and j = A, C, E, while (...) denotes matrix elements which include the proper QCD 
renormalization of the coupling, the masses and the fields. Note that since Q q ± and Qj are 
renormalized operators, all subloop divergences are properly canceled in Eq. (56). Hence 
the renormalization constants Z±j V {p) are sufficient to achieve a finite result. 

The general form of the anomalous dimension tensor for double insertions has been 
derived in [29]. In the following discussion we will only need the explicit expression for the 
part of 7^ v {p) related to anomalous diagrams given by 

ilM = E (^ z kM+(J2^ z kM+J2 iMz^M))^)- ( 58 ) 

k=A,C ^ i=± j=A,C,E 

where and 7fcj(/i) are the elements of the ADM in the (Q%, Q q J) and (Q% Qcs, 

sector, while Z v {p) denotes the renormalization constant of Q v . 

The renormalization constants Zf- v (n) have the following perturbative expansion 

k=l ^ ' 1=0 

for I = P,A. Following the standard MS prescription, Z^-(fi) is given by pure l/e l poles, 
except when % = ± and j = E. In the latter case, the renormalization constant is finite 
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to make sure that the matrix elements of double insertions involving evanescent operators 
vanish in n = 4 dimensions [27-29]. 

The finite parts of Eq. (58) in the limit e going to zero gives the anomalous dimension 
tensor. Performing an expansion in powers of the strong coupling one recognizes that the 
first non-trivial correction to 7^ u (fi) arises at the third order. We obtain 

A(2) _ R 7 A(3,1) , a 7 P(3,1) a 7 (2,1) 7 P(1,0) 9 7 (2,1) y P(l,0) . 7 A(2,0) 7 P(1,1) / Rn x 

which clearly verifies the impact of the finite renormalization of the evanescent operator E q A 
and the axial- vector coupling Q\. Note that in the above equation the superscript A has 
been replaced by P whenever possible. The renormalization constants Z±^'^ and Z±c„ 
are found by calculating the three- and two- loop diagrams shown in Fig. 16, whereas the 
determination of Z± A l)1 J anc ^ ^±Ev requires only a one-loop computation. The relevant 
Feynman graphs are displayed on the left-hand side of Fig. 12. 

On the other hand the pole parts of Eq. (58) must vanish. From this condition one ob- 
tains relations between single, double and triple 1/e poles of the renormalization constants. 
In our case the non-trivial ones read 

(61) 

47 P(3,2) 97 (2,1) 7 P(1,1) _ n v ' 

These equations constitute a powerful check of our three-loop calculation. For instance, 
an erroneous omission of the factor fi~ 2e in the second line of Eqs. (52) changes Eq. (60) as 
well as Eqs. (61), and indeed leads to a failure of the check. 

We now give the values of the renormalization constants entering Eq. (60) for the three 
possible renormalization prescriptions for the open fermion line. The quantities of Eqs. (54) 
and 

Z2ft> = -2 (1 ± C A ) , (62) 

do not depend on the treatment of 75 in the open fermion line. This is not the case for the 
remaining ones. In the MS scheme we find 

Z±A >" ~\-8C F ' Z±C '» ~\3C F ^TC A )^ Z±E: » "to ' (63) 

where the expressions in the first line correspond to the NDR scheme, while the second 
line shows the HV and DRED results. Amazingly, HV and DRED defined as in [51] give 
exactly the same results for all renormalization constants through NNLO. 

Inserting Eqs. (54), (55), (62), and (63) into Eq. (60), we see that 7± i ^' ) vanishes in 
all six renormalization schemes. This non-trivial result implies that anomalous subloops 
involving the Z-boson do not give rise to a NNLO logarithm ln(m^/m|) proportional to 
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rn^/M^, in the decay amplitude of K + — > tx + vv. We have checked the absence of these 
0(a 2 s ) terms explicitly by calculating the three- loop SM diagrams containing an anomalous 
bottom quark loop and verifying that in the limit mj going to zero no IR divergence appear 
in the corresponding amplitude. Beyond NNLO the non- logarithmic pieces of three-loop 
diagrams containing anomalous subgraphs will be relevant and it is highly non-trivial 
whether the cancellation between the effects from top and bottom quark triangles carries 
over to this and higher orders. 



where 



6.5 RG Evolution 

Since in our renormalization scheme specified by the evanescent operators in Eqs. (A. 2) 
the Wilson coefficients C±(fi) evolve independently from each other, Eq. (47) splits into 
two inhomogeneous differential equations. Using Eq. (27) the RG evolution of the Wilson 
coefficients entering the Z-penguin contribution may then be recast into the following 
homogeneous differential equation [9, 58] 

/4C + (/i)\ hM o T £»\ 

dpQi) = 4C_(/i) , 7pM = 7 -(a0 l-M , (65) 

V^W V -rM/ 

which can be solved by the standard techniques [36,46,47] introduced for single operator 
insertions. Since jp(fi) and a s (fi) depend on the number of active quark flavors /, we have 
to solve Eq. (64) separately for fib < \i < fi w and fi c < ji < jib- At the bottom quark 
threshold scale fib additional matching corrections, which will be discussed later in this 
section, have to be taken into account. The Wilson coefficients Cp(/j,) are given by 

Cp(fi) = U p (ii, ii w )C p (ii w ) . (66) 

Keeping the first three terms in the expansions of 7p(/i) and of the QCD /^-function, one 
finds for the evolution matrix Up(fi, fi w ) in the NNLO approximation [39,59] 



Up(fj,,H w ) = kp{^)U { p\^^ w )k P l {ii w ) , (67) 



where 



®s(n) *(i) . fa s (fi)\ 2 ~ (2) 



(68) 
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and 



Up ] 0, £V) = V P diag 



denotes the LO evolution matrix, which is expressed through the eigenvalues a P of ■jp 
and the corresponding diagonalizing matrix Vp: 



(69) 

(0)T 



v P Yp Vp 



(70) 



In order to give the explicit expressions for the matrices J P and Jp we define 



jP = VpS^Vp 1 , 



G 



(fc) 



V P Yp Vp 



(71) 



for k — 1, 2. The entries of the matrix kernels and S*p^ are given by 



-rapSij 



(G«) 



2(3 (l + a p - 



') 



( G g>) 



2ft, (2 + 



•) 



where the first line recalls the familiar NLO result [47], while the second and third represent 
the corresponding NNLO expression derived in [39,59]. 

We will now collect the various expressions that enter the RG analysis of the Z-penguin 
contribution. The LO evolution from fi w down to /x& is described by 



/ 



U 



5(0) 



Vb 


JL J_ 

12 / ^ 23 „ 23 
'lb 







12 
23 



\ 



/ —12 
6/ 23 

13 [Vb 



Vb 



J_ 

, 23 



(73) 



V b J 



Adding an extra index for the number of flavors, the corresponding matrices Jp^ read 



J 



5(1) 



/ 


5165 
3174 





o \ 







2267 







1587 


V 


15857 


15305 


14924 , 
1587 ' 


1587 


3174 



(74) 



35 



and 



J 



5(2) 



/-7.35665 


\ 17.7699 



\ 

-54.9107 
-1.7514 18.3025/ 



(75) 



where in the latter matrix we have employed the numerical value of C(3)- 
The LO evolution from /if, down to /i c is characterized by 



U 



4(0) 



/ 



12 



Vcb 







12 
2.-, 



Vcb 



Va 



_6_ 

11 



Vcb 



Icb 



\ 



Vcb 



where r] cb = a s (/ife)/« s (/i c ). The corresponding matrices Jp take the following form 



(76) 



Vcb J 



J 



4(1) 



/ 3750 


15931 



1875 





3569 
1875 

5427 
1250 



5427 15212 / 

iMfi 1875 ' 






15212 



\ 



(77) 



and 



J 



4(2) 






-50.3422 
-0.3657 






4.91177/ 



(71 



In the last relation terms proportional to ((3) have not been spelled out explicitly again. 
The unbracketed superscripts of the above matrices indicates whether the object belongs 
to the effective theory with five or four active quark flavors. 



6.6 Threshold Corrections 

Since Cp(/x) contains the Wilson coefficients C±(n) it receives a non-trivial 0(a 2 s ) matching 
correction when passing from the effective theory with five active quark flavors to the one 

(2") 

with only four. The explicit expression for the discontinuities 5C± (fib) can be found in 
Eq. (35). In the case of C^(fi) one has to distinguish two possible sources of discontinuities, 
corresponding to the two terms in the first line of Eqs. (45): i) radiative corrections to Q v 
alone and ii) diagrams with double operator insertions (Q±,Qz)- In the first case only a 
single one- and a single two-loop diagram similar to the ones shown on the right of Figs. 7 
and 8 can be drawn. These contributions are canceled by counterterms and the matrix 
elements are zero. Sample diagrams of the second type are displayed on the left of Fig. 12 
and in Fig. 13. Since none of them contains a virtual bottom quark the discontinuities of 
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the corresponding matrix elements vanish identically. By matching the effective theories 
at the bottom quark threshold scale fib we obtain from Eqs. (34) in the NDR scheme 

dC„ [fi b) - 3 ln m 2^ 2645 + 20631 ' V 7935 + 529 Vb ) % 



^1587 6877 
+ f (4 3 ~ Vf 3 ) efGM + ^ - r?f ) tf'w 
where = rrib(mb) denotes the bottom quark MS mass. 
6.7 Matrix Elements 

For scales fi below fi c the transition operator in the case of the Z-penguin contribution is 
simply given by — T p = Ti, p ff with 

K P ff = % 9 a 2ft A c Yl C P (fi)(s Ll ,d L )(u, Ll ^ L ) , (80) 
V2 27rsm z 6> w ' 

which is part of the low-energy effective Hamiltonian that we have encountered already in 
Eq. (1). The bilocal contribution to T p in Eq. (45) has disappeared, because the charm 
quark field is integrated out and the effect from the charm quark loop is absorbed into 
Cp(fi c ) through the matching calculation at fi c . There is still a bilocal contribution from 
the up quark loop, but its contribution is suppressed by a factor of Aq CD /m^ with respect 
to the one stemming from the charm quark. The former is not included in our formalism. 
This power-suppressed contribution contains genuine long-distance effects and has been 
computed in [6]. It will be included in our numerical analysis presented in Sec. 9. 

The local operator entering Eq. (80) has zero anomalous dimension. Therefore Cp(fi) 
is /i-independent for scales ft below ft c . Since we do not need to solve a RG equation 
for fi < fi c , there is no need to express the result in terms of the a s of the effective 
theory containing three active quark flavors and we can avoid to include the non-trivial 
matching corrections of Eq. (32) at the charm quark threshold scale fi c . In terms of the 
a s of the effective theory with four active quark flavors the product Cp(fi c ) (Q u (fi c )) takes 
the following form 

CpW (Q,) (0) = (c^c) (Q,) {0) + ^Cf(fi c ) (QTOio))) , (81) 
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Figure 18: Typical diagrams contributing to the matrix elements of Q± at 0(a s ) (left) and 
0{a 2 s ) (right). 



where we have made use of the fact that the renormalized matrix element of Q v does not 
receive radiative corrections to all orders in a s . Notice furthermore the factor of 4 which 
is a result of our definition of Q± and Q v given in Eqs. (19) and (21). 

In order to complete the evaluation of Cp(fi c ) the renormalized matrix elements of the 
bilocal composite operators Q± are needed. Including corrections up to NNLO, we write 
them in terms of the tree-level matrix element (Q u )^ i n the following way 



Q P M) = ^r^(p c ) + ™ r ™ W m m , (82) 




where r±^\fi c ) and r± (fi c ) codify the one- and two- loop corrections, respectively. Like 
in the case 7± [/ (/i) it is useful to decompose the coefficients further into 

r P J k \, c ) = -~ri (fe Ve) - Q - § sin 2 9^ r V J k \, c ) , (83) 

where the superscript A and V marks the corrections arising from diagrams with a double 
operator insertion (Q±,Q 9 A ) and (Q±,Qy). 

Regulating spurious IR divergences dimensionally we obtain after setting Ca = 3 and 
C F = 4/3 in the NDR scheme 



^ (1 V c ) = -2(l±3)(l-ln4) , 4^(^ = 



m 2 



i {2) (/i c ) = -2 (2 ± 9) - 8 (4 ± 1) In 4 ± 24 In 2 4 , 4 {2) (/i c ) = 



54) 



m 2 c m 2 



where m c = m c (/i c ) denotes the charm quark MS mass. The first line of the above equations 
agrees with the known NLO results [9] if one takes the normalizations of the operators 
and of in Eq. (83) into account. See the discussion after Eq. (51). The second line 
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Figure 19: Typical diagrams contributing to the 0(a 2 s ) matrix element of Q± (left) and 
the 0(a s ) flavor off-diagonal wave function renormalization (right). 

represents the new NNLO expressions. We emphasize that also at the NNLO level only 
the part of the double operator insertion (Q±> Qz) proportional to Q C ± Q C A — Q^±Q\ gives a 
non-zero contribution to the matrix element of Q±. To understand this feature we again 
distinguish Feynman graphs with coupling of Qz to up-type quarks as in Fig. 18 from 
those with down-type coupling of Qz as on the left of Fig. 19. The former diagrams do 
not contribute to the matrix elements at all. In addition the latter diagrams do not arise 
at NLO while at NNLO they give a finite contribution. These terms are again cancelled 
by corrections involving a flavor off-diagonal wave function renormalization as shown on 
the right of Fig. 19. Finally we remark that diagrams with coupling to down-type quarks 
and an insertion of Q\ differ only by a sign from those with an insertion of Q v , because 
Q q ± contains only left-handed down and strange quark fields. We recall that in the full 
theory electromagnetic gauge invariance requires that terms proportional to sin 2 6 W in the 
s — > dvv amplitude add to zero in the limit of vanishing external momenta. The fact that 
the vector part of Q z does not contribute to P C (X) is thus nothing else but the realization 
of this Ward identity in the effective theory. 

6.8 Final Result 

Having calculated all the necessary ingredients of the RG analysis at the NNLO level, we 
are now in a position to present the final result for the Z-penguin contribution to X i (x c ) 
that enters the definition of P C {X) in Eqs. (2) and (114). 

In the NNLO approximation the Wilson coefficient Cp(fi c ) of Eq. (65) has the following 
perturbative expansion 

3rM = 4°> W + « + («4fs>) 2 C<?Vc) , (85) 

where 
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4 1) (^) = ^ (1) f>P (0) f>P (0) 4 0) (^) 

+ v cb u A P {0) (j^-j^u^cPm 

4 2 ) (^ = J^^^^^4°)(^) (86) 

+ Wp (1) ^ (0) - c/p (0) 4°Vw) 

+ Wc6 jJ«^ ^ ((7«(/v) - J P (1) 4 0) (^)) 

+ ^t>p (0) (j P {2) - J P {2) - J P (1) - - 54 2) (^)) ul m cf^ w ) 
+ «# (0) ^ (0) (4 2) (M - ip (1) 4 1} (^) - (^ (2) - (Jp (1) ) 2 ) 



The explicit expressions for Cp\fi w ), Up , Jp , and 5Cp\fib) have been given in 
Eqs. (73) to (79). We also recall that rjb = a s (fj, w ) / ol s {^) and r\ ch = (y s (fib)/<^s{^c)- 

It is useful to express the running charm quark MS mass m c (/x) entering the definition 
of the dimension-eight operator Q v of Eq. (21) in terms of the input parameter m c (m c ). 
At the scale \i c the required NNLO relation reads 

tm> = * (i + + (^) V>) i.w . (87) 

Here k c = 7/c 4//25 with ?7 C = a s (/i c )/a s (m c ) and 

(1) 15212 , _ „ 

^ 5 (88) 

(2) _ 966966391 231404944 _ x 272751559 _ 2 128 , _ 2 . 

c " 10546875 ~ 3515625 Vc " 10546875 Vc Vc ' ^ ' ' 

Including corrections up to third order in perturbation theory one finds from Eq. (81): 

CM = Ko ^ (d£) c!?) W + c? w + w) • < 89 > 

where 

4 0) (^) = ^ P(0) (^), 

4 1} W = OO + 4 £ cf (0) („ c ) pf (1) („ c ) + ) ^ P(0) 00 , 
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4% c ) = C^\p c ) + 4^ {c^\p c )p^\p c ) + C^\p c )p^\p c )) 

i=± 

+ & (V W (Aic) + 4 £ Cf <°> („ c ) pf (1) („ c )) + a P(0) (AO • (90) 



The factors of 4 are again a result of the definition of Q± and Q;, in Eqs. (19) and (21). 

The coefficients p±^ k \p c ) are obtained from r±^ k \p c ) by expanding the charm quark 
MS mass m c (p c ) entering Eqs. (84) in a s around m c (m c ). Explicitly one finds 

p± (1) (Pc) = r± (1) (/i c ,m c ) + (1 ± 3) hxK c , 

a 2 ( 91 ) 
p± (2) (/i c ) = r£ (2) (/i c , m c ) - 4 (4 ± 1) In k c t 12 In 2 Kc ± 24 In k c In ^| + (1 ± 3) . 

m c 

Here the additional argument in r± \p c ,m c ) indicates that the expansion coefficients of 
Eqs. (83) have to be evaluated at m c (m c ) and not at m c (p c ). Note that the second term 
in the first line of the above equations is absent in the analytic NLO formulas of P C {X) 
presented in [9,10]. The relevance of this /independent term will be discussed in Sec. 9. 



7 Electroweak Box Contributions 
7.1 Effective Hamiltonian 

Apart from the presence of a non-trivial matching correction at the bottom quark threshold 
scale pi, the NNLO correction in the electroweak box sector does not involve new conceptual 
features compared to the LO and NLO. This simplifies the following discussion notably. 

For scales p in the range p c < p < p w the electroweak box contribution involves 
defined in Eq. (42) as well as the effective Hamiltonian given by 

H % = % K aJ™r C ?MQ» ■ 02) 
V2 M£ sm z 9 W 

In terms of these building blocks the part of the transition operator T stemming from 
electroweak boxes can be written as 



-T 



B 



(93) 



Notice that in passing from the first to the second line we used the fact that the Wilson 
coefficients of the effective charged-current couplings Ql and Q\ equal one at all scales. In 
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Figure 20: Examples of electroweak box diagrams that contribute to the 0(a 2 s ) correction 
to the initial value of the Wilson coefficient of the leading dimension-eight operator Q v . 

particular they do not receive matching corrections at any scale. The last term in Eqs. (93) 
is the bilocal composite operators Q B which has been introduced in Eq. (20) already The 
factor of 4 originates once again from the use of chiral fermion fields in Q B and Q v . 



7.2 Initial Conditions 

The initial conditions of C B (fi) are as usual found by matching perturbative amplitudes 
in the full and the effective theory. Examples of two-loop electroweak box diagrams can 
be seen in Fig. 20. Regulating spurious IR divergences once dimensionally and once with 
small quark masses we found identical results for the initial conditions. In the NDR scheme 
supplemented by the definition of the evanescent operator given in Eq. (A. 3), the non-zero 
matching conditions read 



C ! B(1) (^ ) = -4(9 + 41n^ , 

C? (a W) = -8CV (20 + 2vr 2 + 25 In || + 6 In 2 ||) 



(94) 



where the first line agrees with the literature [9,10], while the second one is the new NNLO 
expression. 7 



7.3 Anomalous Dimensions 

In the case of the Wilson coefficient C^{p) the RG equation takes the following form 

/^^(/i) = T.MCfW + 4 7i ? (jm) , (95) 

We remark that the logarithmic term in the second line of Eqs. (94) differs from the expression one 
would expect from the results published in [9] . The disagreement is due to a subtlety in regulating spurious 
IR divergences [10,13]. This mistake has been corrected in [10]. 
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Figure 21: Typical examples of diagrams that describe the mixing of the double insertion 
(Qs,QD into the leading dimension-eight operators Q v at 0(ot 2 s ) (left) and 0(a^) (right). 



with 7„()ii) given in Eq. (48). The anomalous dimension tensor 7^(/x) encodes the mixing 
of the bilocal composite structures Q B into Q v . Sample diagrams are shown in Fig. 21. 
The UV pole parts of these Feynman graphs are evaluated using the method that has been 
described earlier. The factor of 4 in the above equation is a direct result of the factor of 4 
in Eq. (93). 

In the NDR scheme supplemented by the definition of the evanescent operator given in 
Eq. (A. 3) the expansion coefficients of J^ifj) read 

7* (0) = "8 

7 i ? (1) = 8CV, (96) 

^ = ^g-fc.-(£-^)c(3) + f/), 

where the second line differs from the findings for 7^ of the original NLO calculation [9] 
even after taking into account a factor of 4 stemming from the different normalization of 
Q B and Q v used here and therein. It however agrees with the results of [10] where the 
error made in [9] has been corrected. The third line represents our new NNLO result. 



7.4 RG Evolution 

Obviously the RG evolution of the electroweak box contribution may be recast into the 
following homogeneous differential equation 

= 1 T b^)Cb^) , (97) 

where 

4 "\ *. ,.A_fO 
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As both 7s(/x) and a s (fi) depend on the number of active quark flavors /, we have to 
solve Eq. (97) separately for [l w and [i c < fi < fi b . At the bottom quark threshold 

scale /if, additional matching corrections, which will be discussed in the next subsection, 
have to be taken into account. 

In the following we will detail the different expressions that enter the RG analysis of the 
electroweak box contribution. Our notation derives from Eqs. (67) to (72) thereby. The 
LO evolution from \x w down to /x& is related to 



U 



5(0) 
B 
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The corresponding matrices Jg are given by 



J 



5(1) 



B 



1 






Vb a ) ' 


by 










2402 


14924 


1587 


1587 



(99) 



and 



jT = I 











1296371522 _ 34624 a/o\ 177621017 , 800 /■ (o\ 
\ 39457581 1081 7555707 ' 23^^% 

The LO evolution from \x\, down to [i c is induced by 

1 
v -12(l -J?"*) ^ 

The corresponding matrices read 



(100) 



(101) 



u 



4(0) 



(102) 



J 



4(1) 









581 


15212 


1875 


1875 



and 



J 



4(2) 
B 



( 



\ 684990354 _ 6976 /-(o\ 
\ 19140625 245 ^ v ' 







272751559 , 128 



10546875 1 5 



+ fC(3), 



(103) 



(104) 



7.5 Threshold Corrections 

In analogy to the case of C^(fi) all discontinuities of the matrix elements that could 
potentially contribute to the threshold correction of C^(n) vanish identically. By matching 
the effective theories at the bottom quark threshold scale we obtain from Eqs. (34) in the 
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Figure 22: Typical examples of diagrams that contribute to the matrix elements of Q at 
0(a s ) (left diagram) and 0{a 2 s ) (right diagram). 



NDR scheme supplemented by the definition of the evanescent operator given in Eq. (A. 3) 
the following non-trivial correction 

^w.-f^K^- W**-™ + * B<?MW ) ■ (105) 

where m& = mb(rrib) is the bottom quark MS mass. 
7.6 Matrix Elements 

For scales \i below \i c the transition operator in the case of the electroweak box contributions 
takes the form — T B = 7i B s with 

Kl = % 9 *« C l M{sa,d L ){y tLl ^ L ) , (106) 

V2 27rsm 9 W t *-^ 

which is part of the low-energy effective Hamiltonian of Eq. (1). Again the bilocal contri- 
bution to T B in Eq. (93) has disappeared, because the charm quark field is integrated out 
and the effect from its loop is absorbed into C^(/i c ). The leftover contribution from the up 
quark loop to T B is like in T p power-suppressed. The numerical size of these corrections 
has been calculated in [6] and is included in our numerical analysis of Sec. 9. 

It is again favorable to express the final low-energy Wilson coefficient in terms of the 
a s of the effective theory with four active quark flavors rather than in terms of the a s of 
the effective theory with three active quark flavors. Proceeding in this way the product 
Cg(^ c ) (Qu(^c)) can be written in the form 

C B M (Q.} i0) = ^jr^T (Q,) (0) + 4 (Qf (*,)>) , (107) 
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where Qf denotes the part of Q B that contains a lepton of flavor i. The factor of 4 is 
again a consequence of the definition of Q B and Q u in Eqs. (20) and (21). 

The renormalized matrix elements (Qf(fi c )) are found by computing the finite parts 
of one- and two-loop diagrams. Examples can be seen in Fig. 22. Regulating spurious 
IR divergences once dimensionally and once with small quark masses we found identical 
results for the matrix elements. In the NDR scheme supplemented by the definition of 
the evanescent operator given in Eq. (A. 3) the expansion coefficients for the tau lepton 
contribution to the matrix element take the following form 

rfC 1 )( Mc ) = 5 + -^ln* T + 41n4, 
1 — x T 

,?» W = -2C F (^-^ + Xt ^^ 3 ) T ) - Y^Li 2 (1 - x T ) (108) 

_n-n^_^24 lni v4_ 61n2 4\. 

V 1 - x T (i-x T ) J m z c m l c ) 

Here x T = vn? T jm\ with m c = m c (/i c ) denotes the ratio of the tau lepton and the charm 
quark MS mass squared. The first line of the above equations agrees with the known NLO 
result [9], after including the aforementioned factor of 4, while the second one represents 
the new NNLO expression. 

In the case of the electron and the muon the lepton mass can be neglected compared 
to the charm quark mass. In the limit x T going to zero Eqs. (108) simplify to 

mj. 

(109) 

^(*)-2C,(9- to g- 6ta 'g). 
7.7 Final Result 

The analytic expression for the electroweak box contribution to X e (x c ) that enters the 
definition of P C {X) in Eq. (2) can be obtained in a straightforward manner following the 
detailed exposition presented at the end of Sec. 6. 

In particular the perturbative expansion of the Wilson coefficient Cs(/i c ) in Eq. (98) is 
given by Eqs. (85) and (86) after replacing all subscripts P by B. The explicit expressions 
for C^ifiw), U f B { °\ J f B (k \ and \fi b ) can be found in Eqs. (99) to (105). 

Keeping terms up to third order in the strong coupling expansion Eq. (107) can be 
written as 

CM = * ^ {-^k) '^ + ^ w + TrPS" w) , mo) 
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where k c is defined after Eq. (87) and 

4 (1 Vc) = C^\^ c ) + Apf W^) + °>(„ c ) , (111) 

4 (2 Vc) = ^ (2) (^) + 4 P f (2) (/i c ) + <7* W^) + e ] C^\» c ) . 

Here the factors of 4 arise again from the use of chiral fermion fields in Q B and Q u of 
Eqs. (20) and (21). 

The coefficients pf^ k \fj, c ) are obtained from rf (/i c ) by expanding the charm quark 
MS mass m c (/i c ) entering Eqs. (108) and (109) in a s around m c (m c ). Explicitly we find in 
the case of the tau lepton 

pf (1) (/ic) = rf (1) (/i c ,™ c ) + — ^— (K c \nK c lnx \ j 

pf (2) (p c )=rf (2) (/i c ,m c )+ 32 ^^"^ 



In x T 



(sc T - k c ) V 3 (1 - x r ) 
x r (x T (13 - 29x r ) + n c (3 + 29x^) - k\ (3 + 13x T )) 
12 (x T - k c ) (1 - x T ) 2 

+ 2x 

H — in k c H In x T In k c — in k c 

b ^-Li2 (1 - x T ) - x T Li 2 ( 1 - — 



2(l-x T ) 



(112) 



32 / k c (1 - x T ) (2x T - K c ) 
+ 7 Typ, r I x T (1 - « c ) H ■ lnx 7 

[X T K c ) [I X T ) \ X T K c 

X(l-«c)(l-2x r + «e) \ /i c 2 



(x T - k c ) (1 - x T ) J m 



+ ^ C ( 1 3— lns T + ^ ln« c Kf. 

In the limit x T going to zero which is the relevant one in the case of the electron and 
the muon one arrives at 



p5 1) (Pc)=r5 1 )(/i c ,m c )-41n 



2 , , (113) 

c 

- L, 1 — — - - L, " - L- fl 

3 777^ 



P5 2) (Pc) = r e B f (/i c , m c ) - I In k c + 16 In 2 k c - 32 In « c In *± - 4£« 



The additional argument in rf^ k \p Cl m c ) signals that the expansion coefficients of 
Eqs. (108) and (109) have to be evaluated at m c (m c ) and not at m c (/j c ). The same applies 
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to the variables x T and m c appearing explicitly in Eqs. (112) and (113). We remark that the 
second terms in the first lines of the latter equations are not present in the NLO formulas 
of P C {X) given in [9, 10]. The importance of this /x c -dependent terms will be discussed in 
Sec. 9. 



8 Final Result for P C (X) at NNLO 

The function X e (x c ) that enters the definition of P C {X) in Eq. (2) is given in terms of the 
contribution of the Z-penguin and the electroweak boxes by 

X £ (x c ) = C P (fx c ) + C e M, (114) 

where the analytic NNLO expression for Cp(/x c ) and Cg(/i c ) can be found in Eqs. (89) and 
(110), respectively. The latter equations together with Eqs. (2) and (114) then enable one 
to find the analytic formula for P C (X) through 0(a s ). 

The explicit analytic expression for P C (X) including the complete NNLO corrections 
is so complicated and long that we derive an approximate formula. Setting A = 0.2248, 
nT'tlrrit) = 163.0 GeV and fi w = 80.0 GeV we derive an approximate formula for P C (X) 
that summarizes the dominant parametric and theoretical uncertainties due to m c (m c ), 
a s (M z ), /i c , and fib- It reads 

/ f \ \ 1-3750 / /A/T \ \ L9480 / \ 
PJX) = 0.3832 ( mc[ y ) ( ) 1 + V K ijlm Ll V a L k L l , (115) 

cV ; Vi- 30GeV / v o.ii87 y \ f-^ l 3 c s Mc Mb / 



where 



L = In I mc{mc) ) , L n = In (^ Mz 



1.30 GeVy 3 V 0.H87 

Mc V 1.5 GeV/ M6 V5.0GeV/ ' 



;il6) 



and the sum includes the expansion coefficients given in Tab. 3. The above for- 
mula approximates the exact NNLO result with an accuracy of ±0.6% in the ranges 
1.15 GeV < m c (m c ) < 1.45 GeV, 0.1150 < a s (M z ) < 0.1230, 1.0 GeV < \i c < 3.0 GeV, 
and 2.5 GeV < /i& < 10.0 GeV. The uncertainties due to m t (m t ), fi w , and the different 
methods of computing a s (/x c ) from a s (M z ), which are not quantified above, are all below 
±0.2%. Their actual size at NLO and NNLO will be discussed in the next section. 



18 



n woo = 0.7432 


^oioo = -3.3790 


^ooio = -0.0001 


Koooi = 0.0028 


Knoo = -0.4350 


Kioio — 0.1669 


K0020 = -0.0903 


«oo2o = -0.0065 


^0030 = 0.03 30 



Table 3: The coefficients Kijki arising in the approximate formula for P C {X) at NNLO. 

9 Numerical Analysis 

9.1 Theoretical Uncertainties of P C (X) at NLO 

Before presenting the numerical analysis of the NNLO correction to P C (X), it is instructive 
to display the theoretical uncertainties present at the NLO level. These originate in the 
leftover unphysical dependence on /i c but are also due to the dependence on fi b and fi w , and 
to higher order terms that arise in the evaluation of a s (fi c ) from the experimental input 
a s (M z ). The latter uncertainties have not been scrutinized in previous NLO analyses of 
the charm quark contribution [9,10]. 

The dependence of P C {X) on fi c can be seen in Fig. 23. The solid red line in the upper 
plot shows the NLO result obtained by evaluating a s (fi c ) from a s (M z ) solving the RG 
equation of a s numerically, while the long- green and short-dashed blue lines are obtained 
by first determining the scale parameter A^jg from a s (M z ), either using the explicit solution 
of the RG equation of a s or by solving the RG equation of a s iteratively for Aj^g, and 
subsequently calculating a s (fj, c ) from Aj^g. The corresponding two-loop values for a s (/x c ) 
have been obtained with the program RunDec [60]. Obviously, the difference between the 
three curves is due to higher order terms and has to be regarded as part of the theoretical 
error. With its size of ±0.012 it is comparable to the variation of the NLO result due to 
fi c , amounting to ±0.020. 

In [9, 10] values for the latter uncertainty have been quoted that are more than twice 
as large. The observed difference is related to the definition of the charm quark mass. 
While in the latter publications the value m c (m c ) has been employed in the logarithms 
ln(jil/ml) of the one-loop matrix elements, we consistently apply m c (/i c ) throughout our 
NLO analysis. 

Using the MS scheme and integrating out the charm quark field at the scale /x c , the mass 
m c (n c ) appears in all intermediate steps of the computation. In particular, the bare one- 
loop matrix elements computed at NLO are proportional to m 2_2e (/i c ). After multiplying 
this with the 1/e poles and finite parts of the loop integrals, adding the counterterm 
diagrams and expanding in e one finds the results given in the first lines of Eqs. (84), (108), 
and (109). Switching by hand to m c (m c ) in the argument of the logarithms ln(/i 2 /m 2 ) 
as done in [9], amounts thus to treating the factors of m 2 (/i c ) and m~ 2e (yU c ) differently, 
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Parameter 


Value ± Error 


T3 „f „„„„„„ 

Keterence 


m c (m c ) [GeV] 
a s (M z ) 
lm\ t [lO" 4 ] 
ReA t [10" 4 ] 
ReA c 


1.30 ±0.05 
0.1187 ±0.0020 

l.tVI _ .098 
q i O+0.20 

-0 22006+ 000093 

u - zzuuu -0. 00091 


[61], our average 
[24] 
[19] 
[19] 
[19] 



Table 4: Input parameters used in the numerical analysis of P C (X), B(K + —> 7r + ^z/), | V*rf 
sin 2/5, and 7. 



although they stem from the very same analytical term in n = 4 — 2e dimensions. It 
is important to realize that such a replacement i) introduces a correction which has no 
diagrammatic counterpart at two loops, and ii) implies a mass counterterm at NNLO 
which is not MS, as it contains an explicit \n(fi 2 /m 2 ) term. Using m c (m c ) in the \Yi(^ 2 /m 2 ) 
terms of the one-loop matrix elements is for these reasons disputable, although it leads to 
results that differ from the one obtained with m c (// c ) by terms that are formally of NNLO. 
These 0(a s ) terms lead to an artificially large /^-dependence at NLO and are hence not 
a good estimate of the size of the uncalculated higher order terms: substituting all factors 
m c (/j, c ) entering the one-loop matrix elements in a consistent way by m c (m c ) results in an 
uncertainty from \i c which is close to the one quoted above. In practice this replacement is 
achieved by employing m c (m c ) in the logarithms \n(n 2 Jm 2 c ) of Eqs. (84), (108), and (109), 
and by rescaling the NLO terms stemming from the second terms in Eqs. (81) and (107) 
by a factor k" 1 . Of course our conclusion is also based on the actual NNLO calculation, 
which indeed finds a correction that is much smaller than the theoretical uncertainty at 
NLO reported in [9,10]. 

Finally, while in [9, 10] only /i c was varied, the theoretical error given in Eq. (3) in- 
cludes also the dependence on /j, b and fi w of ±0.001 each. The specified scale uncertain- 
ties correspond to the ranges 1.0 GeV < /x c < 3.0 GeV, 2.5 GeV < /x& < 10.0 GeV, and 
40.0 GeV < /jl w < 160.0 GeV, and the quoted theoretical error has been obtained by vary- 
ing them independently. 

9.2 Branching Ratio for K + — > ir+vv at NLO 

Using the input parameters listed in Tab. 4, we find from Eqs. (3), (6), (7), and (12) at the 
NLO level 

B(K + ti + vu) = (7.96 ± 0.76 Pc{ x) ± 0.84 othcr ) x 10" 11 , (117) 
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Figure 23: P C {X) as a function of \i c at NLO (upper plot) and NNLO (lower plot). The 
three different lines correspond to three different methods of computing a s (fi c ) from a s (M z ) 
(see text). 
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where the second error collects the uncertainties due to k+, 5P C)U , X(xt), and the CKM 
elements. Numerically, the enhancement of B(K + — > tt + uu) coming from 8P CyU [6] has been 
compensated by the suppression due to the decrease of M t [21]. 

9.3 Branching Ratio for Kl — ► tt q uv at NLO 

Employing the value and errors of ImAj as given in Tab. 4, we obtain from Eqs. (12), (17), 
and (18) in the NLO approximation 

B(K L -> ttVz/) = (2.85 ± 0.05^ ± 0.39 other ) x 10" 11 , (118) 

where the second error collects the uncertainties due to kl, rnt{ m t), an d the CKM elements. 

9.4 Theoretical Uncertainties of P C (X) at NNLO 

Having described the details of our calculation in the previous sections, we now present 
our results for P C (X). From Eqs. (89), (110), and (114) we find at NNLO 

P C (X) = 0.375 ± 0.009 thcory ± 0.031 mc ± 0.009 aa . (119) 

Obviously the error on the charm quark contribution to K + — > tx + vv is now fully dominated 
by the uncertainty in m c (m c ). Comparing these numbers with Eq. (3) we observe that our 
NNLO calculation reduces the theoretical uncertainty by a factor of 4. 

As can be nicely seen in Tab. 1 and in the lower plot of Fig. 23, P C {X) depends very 
weakly on \i c at NNLO, varying by only ±0.006. Furthermore, the three different treat- 
ments of a s affect the NNLO result by as little as ±0.001. The three-loop values of a s (fj, c ) 
used in the numerical analysis have been obtained with the program RunDec [60]. The 
theoretical error quoted in Eq. (119) includes also the dependence on \x h and \x w of ±0.001 
each. The presented scale uncertainties correspond to the ranges given earlier, and the 
specified theoretical error has again been obtained by varying them independently. 

9.5 Branching Ratio for K + — > 7r + ui> at NNLO 

Using Eqs. (6), (7), (12), and (119) the result in Eq. (117) is modified to the NNLO value 
B(K + -> tx+vv) = (8.01 ± 0A9 Pci x) ± 0.84 othcr ) x 10~ n . (120) 

As the value of m c (m c ) is, besides the CKM parameters, the main leftover parametric 
uncertainty in the evaluation of B(K + — > tt + uu), we show in the upper plot of Fig. 3, 
B(K + — > it + uu) as a function of m c (m c ). At present the errors from the CKM parameters 
veils the benefit of the NNLO calculation of P C (X) presented in this paper. For complete- 
ness in the lower plot of the same figure we also show the dependence on M t , which is 
significantly smaller, because M t is much better known than m c (m c ). 
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Observable 


Central ± 68% CL 


95% CL 


99% CL 


B(K L -► 7T°uu) [10- 11 ] 


9 77 +0.58 
z -' '-0.49 


[1.80,3.82] 


[1.62,4.23] 


Pc(X) 


n o fi7 +o.050 

u -°°' -0.049 


[0.306, 0.427] 


[0.293,0.439] 


B{K + -» 7r+i/i>) [10~ n ] 


7QK+1.81 
'• yo -1.67 


[5.61,10.50] 


[5.28,11.18] 




2.56 ±0.30 


[2.02,3.16] 


[1.88,3.38] 




0.375 ± 0.024 


[0.333,0.418] 


[0.324,0.429] 




7.68 ±0.58 


[6.61,8.87] 


[6.32,9.26] 



Table 5: SM predictions for B(K L -> ir vv), P C (X), and B(K + -> it+uv) at the 68%, 
95%, and 99% CL obtained by the CKMfitter Group (upper half) and the XJTfit Collab- 
oration (lower half) incorporating all constraints on the CKM elements following from a 
global analysis of the standard UT. 



9.6 Statistical Analyses of the Branching Ratios of K — ► izvv 

The partial uncertainties given in Eqs. (118) and (120) are not statistically distributed. A 
very important issue in determining the central SM values and errors of B(Kl — > n°uu) and 
B(K + — > 7i + uu) is thus the treatment of the experimental and especially the theoretical 
uncertainties entering these observables. The increasing accuracy in the global analysis of 
the standard UT and the achieved reduction of the theoretical uncertainty of P C (X) clearly 
calls for a closer look at the matter in question. 

To this end it is of interest to see what results are obtained by the two most devel- 
oped statistical methods, namely the iffit approach used by the CKMfitter Group and the 
Bayesian approach employed by the XJTftt Collaboration and to identify those experimental 
and theoretical uncertainties for which a reduction of errors would contribute the most to 
the quality of the determination of the K — > i\vv branching ratios. In this context we 
would like to caution the reader that a direct comparison of the results obtained by the 
two groups in Tab. 5 is quite challenging and the comments given below are hopefully of 
help for the reader to make her or his unbiased judgment of the situation. Our final result 
is given subsequently in Sec. 9.7. 

The numerical results for B(Kl — * ir uu), P C (X), and B(K + — > -k + vv) obtained by the 
CKMfitter Group and the XJTftt Collaboration are summarized in Tab. 5. The corresponding 
likelihood and probability density functions are displayed in Figs. 24 and 25. Apart from 
the CKM elements the employed input agrees with the one that has been used to obtain 
the numerical values for B(Kl tt°ujj), P c (X), and B(K + — > tt + vv) presented earlier in 
Eqs. (118), (119) and (120). 

While the CKMfitter Group and the XJTftt Collaboration find comparable errors at 95% 
and higher confidence levels (CL), the XJTfit Collaboration obtains significantly smaller 
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Figure 24: Likelihood of B(Kl — > 7r uu), P C (X), and B(K + — > ii + vv) following from a 
global CKM analysis performed by the CKMfitter Group. The borders of the 68% (95%) 
probability regions are indicated by the long (short) dashed lines. 
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B(K+ -> tt+vv) [KT 10 ] 



Figure 25: Probability density functions of B{Kl — > ir Q vv), -P c (^0 and — > 7r + z/z/) 
following from a global CKM analysis performed by the Ul/ii Collaboration. Dark (light) 
areas correspond to the 68% (95%) probability regions. 
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errors at the 68% CL. This difference is expected, because both groups treat theoretical 
errors differently: the XJTfit Collaboration assigns a probabilistic meaning to them while the 
CKMfitter Group scans for a best fit value. It can be most easily understood by discussing 
the value of P C {X) at NNLO. Here the dominant error is the parametric uncertainty in 
m c (m c ), which is treated as a theoretical uncertainty by both groups. By assigning a flat 
probability density for m c (m c ) within the error range of 1.25 GeV < m c (m c ) < 1.35 GeV 
only 68% of this measure is used by the XJTfit Collaboration to compute P C (X) at the 68% 
CL. In consequence the error of P C {X) found by the XJTfit Collaboration is smaller by a 
factor of 0.68 compared to the error one would obtain by treating the uncertainty in m c (m c ) 
as the la range of a Gaussian distribution. On the other hand the CKMfitter Group uses 
the whole parameter range of m c (m c ) to compute the error on P C (X) independently of 
the CL. For a flat probability density the measure is proportional to the CL such that 
at the 95% CL nearly the whole range of m c (m c ) is used by the XJTfit Collaboration. 
Correspondingly the difference in the results for the errors of the CKMfitter Group and 
the XJTfit Collaboration decreases strongly at the 95% and higher CL. 

A detailed analysis of the individual sources of uncertainty entering the SM prediction 
of B(Kl — > tt°uu) and B(K + — > tt + uu) using a modified version of the CKMfitter package 
leads to the following picture. In both cases residual scale uncertainties are no longer a 
dominant source of error as they numerically amount to around 9% and 11% of the total 
error only. Hence other intrinsic theoretical errors come to fore. In the case of B(Ki — > 
iftvv) the error associated with tk l is now the main source of theoretical uncertainty 
since the error of these isospin breaking corrections of ±3% translates into around 15% 
of the total uncertainty. In the case of B(K + — > tt + uu) the uncertainties associated with 
the parameters 5P C , U and become prominent. Numerically the total error introduced 
by the long-distance and isospin breaking corrections amounts to about 28% of the final 
uncertainty in B(K + — > ir + uu). This error is thus slightly larger than the error due to 
the charm quark mass which for m c (m c ) = (1.30 ± 0.05) GeV amounts to roughly 20% of 
the final uncertainty. The remaining errors of about 76% and 41% for B(Kl — > 7r°z/z/) and 
B(K + — > 7i + uu), respectively, are due to the uncertainty in the top quark mass, a s (M z ), 
and the CKM elements. The given numbers have been obtained by removing the individual 
errors from the fit in the order the have been mentioned in the text. This study makes 
clear that if one wants to achieve predictions of BiK^ — > tt°vP) and B(K + — > tt + vv) at the 
level of ±5% or below further theoretical improvements concerning the isospin breaking 
corrections, long-distance effects, and the determination of the charm quark mass are 
indispensable. 

9.7 Final Predictions for Kl — > 77°ui> and K + — > 7T + i/i> 

Given the sizeable difference of the 68% CL error intervals obtained by the CKMfitter 
Group and the XJTfit Collaboration we base our final predictions on our own analysis of 
the charm quark contribution and the branching ratios. To determine our final results for 
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B(Kl — > ifivv), P c (X), and B(K + — > tx + vv) we could in principle add the errors given 
in Eqs. (118), (119), and (120) linearly. These estimates might be too conservative. On 
the other hand adding the errors in quadrature would be probably too optimistic since 
the uncertainties are not statistically distributed. Therefore we quote as the final result 
the mean of the values obtained by adding the individual errors once linearly and once in 
quadrature. In the case of — > -n^vv this gives 

B(K L -> n°uu) = (2.8 ± 0.4) x 10~ n , (121) 

while in the case of K + — > ir + vv one has 



P C (X) = 0.38 ±0.04, (122) 

and 

B(K + -> tt+z/P) = (8.0 ± 1.1) x 10~ n . (123) 

The given uncertainties represent the ranges in which we believe that the true values of 
B(Kl — > 7r°z/p), P C (X), and B(K + — > tt + vv) are located with high probability. 



9.8 Impact on the Determination of the CKM Parameters 

As seen in Eq. (4) the accuracy of the determination of |Vfd|, sin 2/5, and of the angle 7 in 
the UT depends sensitively on the error in P C (X). 

The reduction of the theoretical error in P C {X) from ±9.8% down to ±2.4% translates 
into the following uncertainties 

g(IKdl) 
\Vu\ 

<7 (sin 2/3) = { ' ' (124) 



±4.0% , 


NLO, 


±1.0%, 


NNLO, 


±0.024, 


NLO, 


±0.006, 


NNLO, 


±4.7° , 


NLO, 


±1.2°, 


NNLO, 



a (7) = 

implying a very significant improvement of the NNLO over the NLO results. In obtaining 
these numbers we have used sin 2(3 = 0.724 and 7 = 58.6° [19], and included only the 
theoretical errors quoted in Eqs. (3) and (119). 

A comparison of sin 2(3 determined from clean 5-physics observables with sin 2(3 inferred 
from the K — > txvv system offers a precise and highly non-trivial test of the CKM picture. 
Both determinations suffer from small theoretical errors and any discrepancy between 
them would signal non-CKM physics. The impact of future accurate measurements of 
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Figure 26: Standard UT from future measurements of B(Kl — > ir°uu) and B(K + — > tt + vv) 
with an accuracy of ±10%. Dark (light) areas correspond to the 68% (95%) probability 
regions. In the upper (lower) panel the 68% and 95% domains following from the K — > ttuu 
constraint (present global CKM fit) as found by the XJTfit Collaboration are overlaid. The 
lower plot also shows for comparison the present constraint coming from ek- 
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Figure 27: Constraints in the p-fj plane arising from future measurements of B(Ki — > ^vv) 
and B(K + — > ir + uu) with an accuracy of ±10%. The upper (lower) panel shows the 
likelihood obtained assuming no (optimistic) improvement of both the theoretical and 
experimental errors of the input parameters. For comparison the 68% probability region 
following from the present global CKM analysis of the CKMfitter Group is also displayed. 
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K+ _> tt+z/P and # L -> 7r°w leading to i3( J ftT+ -> tt+z/z/) = (8.0 ± 0.8) x 10~ n and 
-> 7T°z/z/) = (3.0 ± 0.3) x 10~ n is illustrated in Figs. 27 and 26. As can be seen in 
Fig. 26 the expected precision of the determination of (p, ff) from the K — > t\vv system 
cannot quite compete with the one from the present global CKM fit performed by the XJTfit 
Collaboration. On the other hand assuming a reduction of the errors in M t , m c (m c ), r^+, 
tk l , and SP CtU by a factor of 3 would put the K — > i\vv system and the global CKM fit 
almost on the same level. The great potential of the K — > -nvv system is clearly visible 
in the lower panel in Fig. 27 which shows the constraint in the p-f\ plane obtained by the 
CKMfitter Group adopting this futuristic scenario. 

Obviously, the future of the determination of the standard UT from the K — > i\vv 
system will depend on the uncertainties in the measured branching ratios, on the value of 
Tn c (m c ) and also on \ V c b\. Further theoretical improvement concerning the isospin breaking 
and long-distance corrections would be desirable in this respect too. A corresponding 
numerical analysis can be found in the updated version of [3], where the NNLO correction 
to -P C (A) presented here, will be soon included. 

While the determination of \V t d\, sin 2/3, and 7 from the K — > txvv system is without 
doubt of interest, with the slow progress in measuring the relevant branching ratios and 
much faster progress in the determination of the angle 7 from B s — > DK system to be 
expected at LHC, the role of the K — > nvv system will shift towards the search for new 
physics rather than the determination of the CKM parameters. 

Indeed, determining the CKM parameters from tree diagrams dominated K- and B- 
decays and thus independently of new physics contributions will allow to find the "true" 
values of the CKM parameters and the so-called reference unitarity triangle [62]. Inserting 
these, hopefully accurate, values in the formulas for the branching ratios presented here, 
will allow to obtain very precise preditions for the SM rates of both decays. A comparison 
with future data on these decays will then give a clear signal of potential new physics 
contributions in a theoretical clean enviroment. Even deviations by 20% from the SM 
expectations could be considered as signals of new physics, whereas it is not possible in 
most other decays in which the theoretical uncertainties are at least 10%. 

10 Summary 

In this paper we have calculated the complete NNLO QCD correction of the charm quark 
contribution to the branching ratio for the rare decay K + — > tt + uu in the SM. As the charm 
quark contribution is essentially unaffected by new physics, our results are also valid in 
basically all extensions of the SM. 

The main result of our paper is summarized in the approximate but very accurate 
analytic expression for the relevant parameter P C {X) as a function of m c (m c ) and a s (M z ) 
presented in Eq. (14). The remaining scale uncertainties and the uncertainty due to higher 
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order terms in the computation of a s (/x c ) from a s (M z ), that is sizeable at the NLO level as 
seen in Fig. 23, have been drastically reduced through our calculation to a level that they 
can basically be ignored for all practical purposes. This can be seen in Tab. 1 and Fig. 23. 
Nevertheless, an approximate formula for P C {X) containing the dominant parametric and 
theoretical errors due to m c (m c ), a s (M z ), /i c , and fib has been given in Eq. (115), which 
should be useful for future phenomenological analysis of the rare decay K + — ► -k + vv. 

The values of P C {X) for different m c (m c ) and a s (M z ) are collected in Tab. 2. As 
a s (M z ) is already known with an accuracy of better than ±2%, the main uncertainty in 
the evaluation of P C {X) resides in the value of m c (m c ). Our nominal value for P C (X) = 
0.38±0.04 used in the NNLO prediction for the branching ratio of K + — > 7r + uu in Eq. (123) 
corresponds to m c (m c ) = (1.30 ±0.05) GeV but the master formulas for P C (X) in Eqs. (14) 
and (115) as well as Tab. 2 allow one to calculate P C {X) and B(K + — > tt + uu) for other 
values of m c (m c ). 

With the improved recent evaluation [6] of the long-distance contributions to the charm 
component, that can be further improved by lattice calculations [7], and hopefully an 
increased accuracy on m c {m c ) and r K + in the future, the theoretical computation of the 
relevant decay rate will reach an exceptional degree of precision, subject mainly to the 
uncertainties in the values of the CKM parameters. As the latter errors will be reduced 
to a large extend in the coming years through the i?-decay experiments a prediction for 
B(K + — > 7i + uu) with an accuracy significantly below ±10% will be possible before the end 
of this decade. Such a precision is unique in the field of FCNC processes. 

On the other hand, accurate measurements of B(K + — > 7r + uu), in particular in conjunc- 
tion with B(Kl — > 7r uu), will provide a very important extraction of the CKM parameters 
that compared with the information from 5-decays will offer a truly unique test of the CKM 
mechanism both in the SM and some of its extensions. The drastic reduction of the theo- 
retical uncertainty in P C {X) achieved by the computation presented in this paper will play 
an important role in these efforts and increases the power of the K — > nuu system in the 
search for new physics, in particular if B{K + — > 7r + uu) will not differ much from the SM 
prediction. 

Acknowledgments 

We would like to thank Matthias Steinhauser for providing us with an updated version 
of MATAD [63] and A. Hocker, J. Ocariz, and M. Pierini for useful communications con- 
cerning the global analysis of the UT. Discussions with G. Buchalla, S. Jager, H. Lacker, 
A. Poschenrieder, F. Schwab, L. Silvestrini, D. Stockinger, S. Uhlig, and A. Weiler are 
acknowledged. Finally we are grateful to the CKMfitter Group and the XJTfit Collabora- 
tion for supplying the numbers and plots shown in Sees. 9.6 and 9.8. This work has been 
supported in part by the Bundesministerium fiir Bildung and Forschung under contract 
05HT4WOA/3 and the Schweizer Nationalfonds. U. H. and U. N. have been supported by 



61 



the DOE under contract DE-AC02-76CH03000. 



(A.1) 



(A.2) 



!2 5 



Appendix 

A.l Evanescent Operators 

The evanescent operators that arise as counterterms for the one-, two- and three-loop 
diagrams with insertions of the current-current operators 

Ql = (s Ll ,T a q L ) (q L YT a d L ) , 
Ql = (slJ^l) {q^dL) , 

can be chosen to be 

E\ = (s L l» w jy L )(q L l^T a d L ) - (16 - 4e - 4e 2 ) Q\ , 

El = {sLl^^WLl^^dr) - (16 - 4e - 4e 2 ) Ql , 

/ ^719 \ 

El = (^7M 1 M 2 M3M4M 6 ^°^)(^7 MlM2M3/M/i6 7 1 °^) - [256 - 2246 - —e 2 j Ql , 

El = {s Ll ^ 2m q L ) (q L r m ^ 5 d L ) - (256 - 2246 - ^e 2 ) Q\ 

E\ = (^7M 1 A t2 M3M4M5/,6/, 7 ^ a ^)(gL7 WM2/i3M4At5WM7 T a rf L ) - (4096 - 7680 e ) Q\ 

El = (^7 Atl M 2 M3A t 4M5/,6/, 7 ^)(gL7 MlA12M3M4M5WW ^) - (4096 - 7680e) Q\ . 
Here the shorthand notation 7 M1 ... Mm = 7/ii"""7u m and ■y 1 "^™ = . . . y^™ Deen 
used. The above operators have been defined such that %) the ADM of the operators Q q ± 
introduced in Eq. (23) is diagonal through NNLO, and that ii) their particular structure 
differs only by multiples of e 2 times physical operators from the evanescent operators of the 
"traditional" basis [27,36]. The latter operators can be found by the procedure outlined 
in [64]. Up to three loops they have been given in [39]. Of course, the above choice is not 
unique in the sense, that there are many schemes that would satisfy i) through NNLO. It 
is however not possible to define the set of evanescent operators to be invariant under the 
interchange of color structures and to achieve i) simultaneously. Finally, we remark that 
the e 2 terms of the above one-loop evanescent operators E\ and E\ are unambiguously 
determined by condition i). 

In the case of the electroweak box contribution only one evanescent operator arises 
as a counterterm for the one-, two-, and three-loop diagrams considered in this paper. 
Following [9, 10, 13] we chose it as 

2 

E » = -irkY, (^ M d L )(ve L Y M vu) ~ (16 - 46) Q v . (A.3) 
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Here m c is the charm quark MS mass m c (/x) and the explicit factors g 2 and \x 2e follow 
from the normalization of Q v in Eq. (21). 



A. 2 Change from the "Standard" Basis of Current-Current Op- 
erators 

Beyond LO the anomalous dimensions and the Wilson coefficients depend on the definition 
of the operators in n = 4 — 2e dimensions. So far all of the direct NNLO calculations have 
been performed in the operator basis introduced in [38] , which we will call "standard" basis 
from now on. It consists of the following set of physical and evanescent operators 

Q' T = (Qi, Q'l) , E' T = (E* E'l E'l E'i) , (A.4) 



where 



Qi = (s Llfl T a q L ) (q L rT a d L ) 



L 



E[ q = (s L7m T a qL)(qLl m T a d L ) - 16Q'? 

E' 2 q = (s L7m q L )(q L Y lfl2f * 3 d L ) - lQQ' 2 q , 

E's = (^7 w , 2WWti5 T a 9i )( 9i 7' IlM2 ' iW4M5 T a d i ) - 256Q? - 20E? 

E'Z = {s L l^ 2m5 q L ){qLl^^d L ) - 256Q? - 20i^ , 



(A.5) 



while 



Q T = Q-) , E T = (El El El El) , (A.6) 

denotes the physical and evanescent operators used in this paper. The operators E\ and 
.Eg and their primed counterparts have been omitted in the above equations, because they 
do not affect the change of scheme up to the order considered here. 

In this appendix we demonstrate how the results for the two-loop initial conditions [32] 
and the three-loop ADM [39] of the current-current operators can be transformed to our 
basis. This will serve as a cross-check of our results for C^\jjL w ) and 7± as given in 
Eqs. (25) and (29). In all formulas presented below we set Ca = 3 and Cp = 4/3, while 
the number of active quark flavors / is kept arbitrary. 

The transformations relating the primed with the unprimed sets take the simple form 

Q = RQ>, E = M(E' + e UQ' + e 2 V(f \ , (A.7) 

where the matrix R (M) describes a rotation of the physical (evanescent) operators, while 
the matrix U (V) parameterizes a change of basis that consists of adding multiples of e 
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(e 2 ) times physical operators to the evanescent ones. The matrices introduced in Eq. (A. 7) 
are given by 



R 




M 



1 


20 
. 




1 




\ 



1 



20 1 



/ 



u 
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i 
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4 





\ 







4 











4 












v = 




3712 
25 








144 














V 





144 


) 




V 





8032 
75 


/ 



(A.l 



The change of basis in Eq. (A. 7) is /^-independent and leaves, apart from a global 
rotation, the anomalous dimensions and the Wilson coefficients invariant. It, however, 
induces a finite renormalization of the physical operators. In order to restore the standard 
MS renormalization conditions these contributions must be removed by a change of scheme. 
Hence a e-dependent linear transformation of the operator basis is equivalent to a global 
rotation and a change of scheme [39]. 

The finite renormalization corresponding to the above change of basis can be derived 
with simple algebra. Through 0(a 2 s ) we find 



where 



J QQ 
7(2,0) 



>'(2,2) 
J QE 



-R 



2 QE 



'.(y)yy(o) 



R 



-1 



A/(l,l) A/(l,l) 
^QE ^EE 



1 .7(0)7/(1,1) n 7/(1,1) 



(A.9) 



(A.10) 



As always, the matrix M encoding the rotation of evanescent operators does not affect 
the residual finite renormalization. Eqs. (A.9) and (A. 10) agree with the NLO formulas 
of [29,38] and generalize the NNLO formulas of [39]. 

The only feature which has not been discussed in the literature before is the appear- 
ance of the e 2 -dependent transformation characterized by V in Eq. (A. 7). It induces a 
finite O{of) renormalization described by the second and third term in the second line of 
Eqs. (A.9). The former term is analogous to the first line of Eqs. (A.9). The latter term 
stems from a e-dependent 0(a s ) renormalization. While the e-dependent change of scheme 
removes the finite terms, it still leaves the anomalous dimensions and the Wilson coeffi- 
cients invariant. Yet it induces a finite 0(a 2 s ) renormalization. The transformation to the 
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standard MS scheme definitions requires the removal of this finite term, which is achieved 
by the last term in the second line of Eqs. (A. 9). 

The matrices entering the above equations are found from one- and two-loop matrix 
elements of physical and evanescent operators in the "standard" basis. The 0(a s ) ADM 
takes the following form [38,39,42] 



7 



/(o) 



-4 s 

^ 3 

12 



(A.n; 



For the 0(a s ) mixing into the evanescent operators one obtains [38,39,42] 



J QE 



5_ 
12 




J EE 
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-7 


4 
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-> 

12 


2 
9 






-6 





1 















13 
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28 
9 




I 








-14 


64 
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) 



(A.12) 

It is given 
_^ f (A.13) 

18 J 9 

The general NNLO formulas relating the initial conditions of the Wilson coefficients 
and the ADM in two different schemes have been derived recently [39]. They read 

' ( ° x T 
a s {/j, w 



At 0(a 2 s ) only the mixing of physical into evanescent operators is needed, 
by [39,42] 

1531 5_f _J 1_ f J_ 35 

288 216-' 72 81-' 384 864 



7(2,1) 
^QE 



119 1_ f 

16 18 J 



Air 



r(l,0) 



47T 



(^- 1 ) T C / (// w ) 



(A.14) 



and 



7 (°) 


= Rj'^R- 1 






7 (1) 








7 (2) 


= R^R- 1 




>(2.°) A(0)' 






2/3 1 2§§ ) + 2/3 



9/ o 7 (1,0) 

QQ '7 

r(l,0)\2 



(A.15) 



v(i.o) a(o) 

QQ '7 



r(l,0) 



Eqs. (A. 8) to (A.15) allow one to transform the results for the initial conditions of the 
Wilson coefficients of the current-current operators and their ADM in the "standard" 
basis to the basis used in this paper. This enables us to verify that the well-established 
NLO results [38,39,42] coincide, after the change of scheme, with the expressions presented 
in the second line of Eqs. (25) and (29). For what concerns the NNLO we confirm both 
the two-loop initial conditions of the Wilson coefficients [32] as well as the corresponding 
three-loop ADM [39]. 
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